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Abstract 



In this thesis we study the dark matter problem with particular reference to a 
candidate particle within the Standard Model: the H dibaryon. We consider as well 
a scenario which aims to connect the dark matter origin to the Baryon Asymmetry 
of the Universe, studying the examples of H and of a Beyond-the-Standard- Model 
particle X. 

Strongly attractive color forces in the flavor singlet channel may lead to a 
tightly bound and compact H dibaryon. We find that the observation of A decays 
from doubly-strange hypernuclei puts a constraint on the H wavefunction which 
is plausibly satisfied. In this case the H is very long-lived as we calculate and an 
absolutely stable H is not excluded. We also show that an H or another compact, 
flavor singlet hadron is unlikely to bind to nuclei, so that experimental bounds 
on exotic isotopes do not exclude their existence. Remarkably, the H appears to 
evade other experimental constraints as well, when account is taken of its expected 
compact spatial wavefunction. 

In order to check whether the i7 is a viable DM candidate, we consider experi- 
ments sensitive to light particles. Taking into account the dark matter interaction 
in the crust above underground detectors we find a window in the exclusion limits 
in the micro-barn, m^2.4 GeV, range. Remarkably, this coincides with the range 
expected for the tightly bound H. Having these constraints in mind wc conclude 
that the is a good DM candidate, but its production with sufficient abundance 
in the Early Universe is challenging. 

Finally, we present a scenario in which dark matter carries (anti-) baryon num- 
ber Bx and which offers a mechanism to generate the baryon asymmetry observed 
in the Universe. If a^^^^ < cr^f'"'^, the X's freeze out at a higher temperature and 
have a larger relic density than X's. If mx ^ 4.5 Bx GeV and the annihilation 

iv 



cross sections differ hj O{10%) or more, this type of scenario naturally explains 
the observed ^^dm — 5 fib- Two examples are given, one involving the H and the 
other invoking an hypothetical beyond the Standard Model candidate X. 
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Chapter 1 

Dark Matter and the Baryon 
Asymmetry of the Universe 

The existence of Dark Matter (DM) is, today, well established. The name "Dark" 
derives from the fact that it is non-luminous and non-absorbing; it can be detected 
(so far) only through its gravitational interaction with ordinary matter. One of 
the first evidence for the presence of DM was the 1933 measurement by F. Zwicky 
of the velocities of galaxies which are part of the gravitationally bound COMA 
cluster, [1]. Zwicky found that galaxies are moving much faster than one would 
expect if they only felt the gravitational attraction from visible nearby objects. 
Nevertheless, the existence of dark matter was not firmly established until the 
1970's when the measurement of the rotational velocity of stars and gas orbiting 
at a distance r from the galactic center was performed. The velocity at a distance r 
scales as V ~ ^/M{r)/r, where M(r) is mass enclosed by the orbit. If measurement 
is performed outside the visible part of galaxy, one would expect M(r) ~ const., 
or ^; ~ ^l\fr. Instead, observations show that v ~ const., implying the existence 
of a dark mass, with radial dependence M~r or p~l/r^, assuming spherical 



1 



symmetry. The existence of dark mass is probed on different scales: velocities 
of stars, gas clouds, globular clusters, or, as we have seen, entire galaxies, are 
larger than one would predict based on the gravitational potential inferred from 
the observed, luminous mass. 

More recent methods for direct detection of DM include measurements of X- 
ray temperature of the hot gas in galaxy clusters, which is proportional to the 
gravitational potential field and therefore to the mass of the cluster, and observa- 
tions of the gravitational lensing of galaxies caused by the mass of a cluster in the 
foreground. 

On a more theoretical basis, the presence of dark matter allows to relate the 
anisotropics of the cosmic microwave background (CMB) and the structures ob- 
served in galaxy and Lyman-a surveys to a common primordial origin in the frame- 
work of the inflationary model. The currently most accurate determination of DM 
and baryonic energy densities comes from global flts of cosmological parameters to 
these observations. For instance, using measurements of the CMB and the spatial 
distribution of galaxies at large scales. [2], one gets the following constraints on the 
ratio of measured energy density to the critical energy density per — SHq/SttGn, 

^DM — PDm/ Per and Qb = Pb/ Per'- 

QoMh'^ = 0.1222 ± 0.009, 

Qbh^ = 0.0232 ± 0.0013, (1.1) 

where h — 0.73 ± 0.03 is the Hubble constant in units of 100 km s~^Mpc~^, 
Hq — h 100 km s~^ Mpc~^. These two numbers are surprisingly similar, 

QomM = 5.27 ±0.49, (1.2) 

even though in conventional theories there is no reason for them to be so close: 
they could differ by many orders of magnitude. In our work we will explore the idea 
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that these two numbers might originate from the same physical process, thereby 
making the value of their ratio natural. 

The nature of the dark matter particle is an unsolved problem. Candidates for 
dark matter must satisfy several, basic, conditions: they should have very weak 
interaction with the electromagnetic radiation, they must have a relic density that 
corresponds to the observed DM energy density, eq. (1.1), and they should have a 
lifetime longer then the hfetime of the universe. They should also be neutral parti- 
cles. The models of structure formation based on the inflation scenario prefer the 
so called "cold" dark matter, i.e. dark matter particles which are non-relativistic 
at the time of galaxy formation. In these models dark matter fluctuations, caused 
by primordial fluctuations in the inflaton field, are responsible for growth of the 
structures observed today. Baryons follow DM fluctuations falling into their grav- 
itational wells where they form astrophysical objects. A relativistic dark matter 
would have a large free- streaming length, below which no structure could form 
while cold dark matter has a free-streaming length small enough to be irrelevant 
for structure formation. 

A particle candidate for DM is beheved not to be provided by the Standard 
Model (SM). Neutrinos have long been believed to be good SM candidates but 
based on recently derived constraints on the neutrino mass, it has been realized that 
they cannot provide enough energy density to be the only dark matter component 
(the current limit is ^ 0.07). Therefore the DM candidate is usually looked 
for in physics beyond the Standard Model: the most important candidates today, 
which satisfy the above requirements and are well motivated from particle physics 
considerations are axions and the lightest supersymmetric particle. 

Axions are pseudo Nambu-Goldstone bosons associated with the spontaneous 
breaking of a Peccei-Quinn, [3,4], U{1) symmetry at scale Ja, introduced to solve 
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the strong CP problem. In general, their mass is inversely proportional to /a, 
as niA = 0.62 10~^ eV (10^° GeV//^). The allowed range for the axion mass 
10~^ ^niA^ 10~^ eV, see for instance [5]. The lower bound derives from the con- 
dition that the energy density in axionic DM is not higher than the observed DM 
density, and the upper bound derives most restrictively from the measurement 
of the super nova neutrino signal (the duration of the SN 1987A neutrino sig- 
nal of a few seconds points to the fact that the new born star cooled mostly by 
neutrinos rather than through an invisible channel, such as axions). The axions 
are produced with the DM abundance only for large values of the decay constant 
fA, which implies that they did not come into thermal equilibrium in the early 
universe. They were produced non-thermally, for example through a vacuum mis- 
alignment mechanism, see [6]. Experimental searches for axionic DM have been 
performed dominantly through the axion to photon coupling. The Lagrangian is 

— * — * 

= QA-y E • BcpA, where (pA is the axion field, and qa-j is coupling whose strength 
is an important parameter in axion models; it permits the conversion of an axion 
into a single real photon in an external electromagnetic field, i.e. a PrimakofT 
interaction. Halo axions may be detected in a microwave cavity experiments by 
their resonant conversion into a quasi-monochromatic microwave signal in a cavity 
permeated by a strong magnetic field. The cavity "Q factor" enhances the conver- 
sion rate on resonance. Currently two experiments searching for axionic DM are 
taking data: one at LLNL in California, [7], and the CARRACK experiment, [8], in 
Kyoto, Japan. Preliminary results of the CARRACK I experiment exclude axions 
with mass in a narrow range around lO/ieV as major component of the galactic 
dark halo for some plausible range of qaj values. This experiment is being up- 
graded to CARRACK II, which intends to probe the range between 2 and 50 /xeV 
with sensitivity to all plausible axion models, if axions form most of DM. 

The lightest supersymmetric particle belongs to the general class of weak inter- 
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acting massive particles, the so called WIMPs. WIMPs are particles with mass of 
the order of 10 GeV to TeV and with cross sections of weak interaction strength, 
(7 ~ Gp m^jj^p. Supersymmetry (SUSY) allows the existence of a R-parity sym- 
metry which implies that the lightest SUSY particle is absolutely stable, offering 
a good candidate for DM. In most SUSY models the hghtest supcrsymmctric par- 
ticle is the neutralino, a linear combination of photino, wino and two higgsinos. 
Under the assumption that WIMPs were in thermal equilibrium after inflation, 
one can calculate the temperature at which their interaction rate with Standard 
Model particles becomes lower than the expansion rate of the universe. At that 
point they decouple from thermal bath and their number in co-moving volume 
stays constant at later times (thcj freeze- out). Freeze-out happens at temperature 
T ~ m/20, almost independent of the particle properties, which means that par- 
ticles are non-relativistic at the time of decouphng, making WIMPs a cold DM 
candidates. The abundance of DM in this scenario is 0,dm — 0.1 ph/{av), which 
surprisingly corresponds to the measured DM density for weak cross section values. 
The fact that their mass/cross section value is in the correct range, together with 
good motivation from SUSY, makes them attractive candidates. 

Today, a large part of the parameter space expected for neutralino has been 
explored: a region of mass 10 GeV ^ m < TeV, from cr ~ 10 mb to cr ~ 10~^pb, 
and a new generation of experiments reaching a ~ 10~^'^pb is planned for the near 
future. Since direct and indirect detection searches have been unsuccessful, the 
nature of dark matter is still unresolved, although not all of the possible SUSY 
and axion parameter range has been explored, see for instance [5,9, 10]. We will 
examine alternative candidates for DM in the thesis, in particular candidates con- 
nected to the Baryon Asymmetry in the Universe, described in detail below. 

Another important open question in our understanding of the Universe today 
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is the observed Baryon Asymmetry of the Universe. The Dirac equation places 
anti-matter places on an equal footing with matter: in the Big Bang cosmological 
model the early epoch should contain a fully mixed state of matter and anti- 
matter. As the Universe expanded and cooled this situation would result in the 
annihilation of matter and anti-matter, leaving a baryon mass density today of 
Ofe ~ 4 10~^^ which is much lower than the observed value Qi, ~ 0.04, eq. (1.1) (or 
a baryon to photon ratio Ub/n^ ~ 10~^*, eight orders of magnitude smaller than 
measured). Evidently some mechanism had to act to introduce the asymmetry in 
the baryon and antibaryon abundances and prevent their complete annihilation. 
The apparent creation of matter in excess of antimatter in the early universe 
is called Baryogcncsis, for a review see, for instance, [11]. It was A. Sakharov 
who first suggested in 1967, [12], that the baryon density might not come from 
initial conditions, but can be understandable in terms of physical processes. He 
enumerated three necessary conditions for baryogenesis: 

1. Baryon number violation: If baryon number (B) is conserved in all reactions, 
then a baryon excess at present can only reflect asymmetric initial conditions. 

2. C and CP violation: Even in the presence of B-violating reactions, if CP is 
conserved, every reaction which produces a particle will be accompanied by 
a reaction which produces its antiparticle at the same rate, so no net baryon 
number could be generated. 

3. Departure from thermal equilibrium: the CPT theorem guarantees equal 
masses for particle and its antiparticle, so in thermal equilibrium the densi- 
ties of particles and antiparticles are equal and again no baryon asymmetry 
could exist. 

Several mechanisms have been proposed to understand the baryon asymmetry. I 
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comment on a few of the most important models below. 

Grand Unified Theory (GUT) scale Baryogenesis ( [13-15]): Grand Unified 
Theories unify the gauge interactions of the strong, weak and electromagnetic in- 
teractions in a single gauge group. The GUT scale is typically of the order 10^^ 
GcV so baryogenesis in this model occurs very early in the history of Universe. 
GUT's generically have baryon-violating reactions, such as proton decay (not yet 
observed) and they have heavy particles whose decays can provide a departure 
from equilibrium. The main objections to this possibility come from inflation and 
reheating models, - the temperature of the universe after reheating in most models 
is well below Mqut- Furthermore, if baryogenesis occurs before inflation it would 
be washed out in the exponential expansion. 

Electroweak Baryogenesis ( [16]): In this model baryogenesis occurs in the SM 
at the Electroweak Phase Transition (EWPT) - this is the era when the Higgs first 
acquired a vacuum expectation value (VEV) and the SM particles acquired masses 
through their interaction with the Higgs field. This transition happened around 
100 GeV. The Standard Model satisfies all of the Sakharov conditions: 

(i) In the SM there are no dimension 4 operators consistent with gauge symme- 
try which violate baryon (B) or lepton number (L). The leading operators 
which violate B are dimension 6 operators, which are suppressed by 0(1/M^) 
and the operators which violate L are dimension 5, suppressed by 0(1/M), 
where M is a scale of high energy physics which violates B or L. Inside the 
SM, B and L currents are not exactly conserved due to the fact that they are 
anomalous. However, the difference between — j'^ is anomaly-free and is an 
exactly conserved quantity in the SM. In perturbation theory these effects go 
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to zero, but in non-abelian gauge theories there are non-perturbative config- 
urations which contribute to the non-conservation of currents. The vacuum 
structure of a Yang-Mills theory has an infinite set of states. In tunneling 
between these states, because of the anomaly, the baryon and lepton numbers 
will change. At zero temperatures tunneling effects are exponentially sup- 
pressed by exp(— 27r/Q;). At finite temperatures this rate should be larger. To 
estimate this rate one can look for the field configuration which corresponds 
to sitting on top of the barrier, a solution of the static equations of motion 
with finite energy, known as a sphaleron. The rate for thermal fiuctuations 
to cross the barrier per unit time and volume should be proportional to the 
Boltzmann factor for this configuration, [16-18], F = T^e~'^^^^^^'^ . At high 
temperature vanishes and the transition gets the form F = a^T'^. 

(ii) CP-violation has been experimentally observed in kaon decays and is present 
in the SM. However, SM CP violation must involve all three generations. 
The lowest order diagram that involves three generations and contributes to 
CP violating processes relevant to baryogenesis is suppressed by 12 Yukawa 
couplings. The CKM CP violation contributes a factor of 10~^° to the amount 
of baryon asymmetry that could arise in the SM and a Beyond the Standard 
Model CP violation is usually invoked. 

(iii) Thermal nonequilibrium is achieved during first-order phase transitions in 
the cooling early universe. In the electroweak theory, there is a transition to 
a phase with massless gauge bosons. It turns out that, for a sufficiently light 
Higgs, this transition is of the first order. A first order transition is not, in 
general, an adiabatic process. As we lower the temperature, the transition 
proceeds by the formation of bubbles. The moving bubble walls are regions 
where the Higgs fields are changing and all of Sakharov's conditions are 
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satisfied. It has been shown that various non-equihbrium processes near the 
wall can produce baryon and lepton numbers, [19,20]. Avoiding the washing 
out of the asymmetry requires that after the phase transition, the sphaleron 
rate should be small compared to the expansion rate of the universe, or, as 
we have seen above, that Mw be large compared to the temperature. This, 
in turn, means that the Higgs expectation value must be large immediately 
after the transition. It turns out that the current lower limit on the Higgs 
boson mass rules out any possibility of a large enough Higgs expectation 
value after the phase transition, at least in the minimal model with a single 
Higgs doublet. 

Any baryon asymmetry produced in the Standard Model is far too small to 
account for observations, the main obstacle being the heaviness of the Higgs, and 
one has to turn to extensions of the Standard Model in order to explain the ob- 
served asymmetry. 

Leptogenesis ( [21]): In the last few years the evidence for neutrino masses has 
become more and more compelling. The most economical way to explain these 
facts is that neutrinos have Majorana masses arising from lepton- number violating 
dimension five operators (permitted if fermion carries no conserved charges) . These 
interactions have the form C — j^LHLH. For M — Mpi the neutrino mass would 
be too small to account for the observed values. The see-saw mechanism provides 
a simple picture of how the lower scale might arise. It assumes that in addition to 
the SM neutrinos, there arc some SM singlet, heavy neutrinos, N. These neutrinos 
could couple to the left handed doublets ul providing the correct mass for the light 
neutrinos. 

What is relevant is that heavy neutrinos N can decay, for example, to both 
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H -\- v and H + jy, breaking the lepton number. CP violation can enter through 
phases in the Yukawa couphngs and mass matrices of the A^'s. At tree-level these 
phases will cancel out, so it is necessary to consider one loop diagrams and to look 
at quantum corrections in which dynamical phases can appear in the amplitudes. 
These decays then produce a net lepton number, and hence a net B—L. The result- 
ing lepton number will be further processed by sphaleron interactions, yielding a 
net lepton and baryon number. Reasonable values of the neutrino parameters give 
asymmetries of the order we seek to explain. However, all parameters needed for 
precise calculations are not measured yet (in the case of ul masses and CP violating 
couplings) and one needs some additional information about the masses of the A'^'s. 

Affteck-Dine Baryogenesis( [22]) In supersymmetric theories, the ordinary quarks 
and leptons are accompanied by scalar fields, which carry baryon and lepton num- 
ber. A coherent field, i.e. a large classical value of such a field, can in principle 
carry a large amount of baryon number. Through interactions with the infiaton 
field CP-violating and B-violating effects can be introduced. As the scalar parti- 
cles decay to fermions, the net baryon number the scalars carry can be converted 
into an ordinary baryon excess. The Affleck-Dine mechanism is also a mechanism 
for dark matter creation. Fluctuations in the scalar quark fields ( "Q-balls" ) are a 
dark matter candidate if they are stable. If they are unstable, they can still decay 
into dark matter. Since the Affleck-Dine mechanism describes the production of 
baryons and dark matter it could provide an explanation of the ratio between dm 
and from first principles. If supersymetry is discovered, given the success of 
inflation theory, the Affleck-Dine scenario will appear quite plausible. 

To summarize, the abundance of baryons and dark matter in our Universe poses 
several challenging puzzles: 
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(i) Why is there a non-zero net nucleon density and what determines its value? 

(ii) What does dark matter consist of? Can it be explained as a SM particle? 

(iii) Is it an accident that the dark matter density is roughly comparable to the 
nucleon density, poM — 5 Pat? 

In the next Chapter we outline the scenario which aims to connect and answer 
the questions above. In Chapters 3, 4 and 5 we focus on the two concrete DM 
candidates in this scenario, one being a particle within the Standard Model, and 
the other is BSM candidate. We comment on their particle physics properties and 
experimental constraints. As we will see, SM candidate is ruled out while the 
Beyond the Standard Model candidate is safe by many orders of magnitude. 
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Chapter 2 
Proposed Scenario 



In most approaches the origins of DM and the BAU are completely unrelated - 
baryon number density ub is obtained from baryogensis models while the number 
density of DM udm derives from relic freeze-out calculations and their values natu- 
rally could differ by many orders of magnitude (sec [23-26] for other papers where 
these two problems arc related). In this Section we propose a new type of scenario, 
in which the observed baryon asymmetry is due to the separation of baryon num- 
ber between ordinary matter and dark matter and not to a net change in the total 
baryon number since the Big Bang, [27]. Thus the abundances of nucleons and 
dark matter are related. The first Sakharov condition is not required, while the 
last two remain essential. We give explicit examples in which anti-baryon number 
is sequestered at temperatures of order 100 MeV. 

The CPT theorem requires that the total interaction rate of any ensemble 
of particles and antiparticles is the same as for the conjugate state in which each 
particle is replaced by its antiparticle and all spins are reversed. However individual 
channels need not have the same rate so, when CP is violated, the annihilation 
rates of the CP reversed systems are not in general equal. A difference in the 
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annihilation cross section, cr^^"''^ < aj^"^^^, means that the freeze out temperature 
for X's {Tx) is lower than for X's (T^). After the X's freeze out, the X's continue 
to annihilate until the temperature drops to Tx, removing Bx antinucleons for 
each X which annihilates. 

Assuming there are no other significant contributions to the DM density, the 
present values Uon, n-ox and rigx are determined in terms of mx, Bx and the 
observables and ^ = 7710 10"^° or p^rit- From WMAP, 
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Vio = G.5t'o% 

^ = 5.27 ±0.49. (2.1) 
Ub 

Given the values of these observables, we can "reverse engineer" the process of 
baryon-number-segregation. 

For brevity, suppose there is only one significant species of DM particle. Let 
us define e = Then the total energy density in X's and X's is 

Pdm = mxnx{l + e). (2.2) 

By hypothesis, the baryon number density in nucleons equals the antibaryon num- 
ber density in X and X's, so 

Bxfixil - e) = (njv - n^) = — ■ (2.3) 

rriN 

Thus 

^DM /l + e\ mx 



(2.4) 

Qfe \ 1 — ey mj^Bx 
As long as the DM particle mass is of the order of hadronic masses and e is not 

too close to 1, this type of scenario naturally accounts for the fact that the DM 

and ordinary matter densities are of the same order of magnitude. Furthermore, 

since > 1 , the DM density in this scenario must be greater than the nucleonic 

density, unless mx < m^Bx, as observed. 
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Given the parameters of our Universe, we can instead write (2.4) as an equation 
for the DM mass 

For low baryon number, Bx = 1 (2), this imphes 

mx<4.5(9)GeV. (2.6) 

If dark matter has other components in addition to the X and X, the X must 
be hghter still. The observed BAU can be due to baryon number sequestration 
with heavy DM only if Bx is very large, e.g., strangelets or Q-balls. However 
segregating the baryon number in such cases is challenging. 

As an existence proof and to focus discussion of the issues, we present two 
concrete scenarios. In the first, X is a particle already postulated in QCD, the 
H dibaryon (uuddss). New particle physics is necessary, however, because the 
CP violation of the Standard Model via the CKM matrix cannot produce the 
required (9(20%) difference in annihilation cross sections, since only the first two 
generations of quarks arc involved. The second scenario postulates a new particle, 
we call X4, with mass ^4.5 GeV, which couples to quarks through dimension- 
6 operators coming from beyond-the-standard-model physics. In this case CP 
violation is naturally large enough, (9(10%), because all three quark generations 
are involved and, in addition, the new interactions in general violate CP. Review 
of particle properties of these candidates is given in Sections 3.2 and 5.1. After 
deducing the properties required of these particles by cosmology we discuss indirect 
(Sections 4.2 and 5.2) and direct (Section 4.1) searches. As we shall show, the H, 
H scenario can already be ruled out by limits on the heat production in Uranus, 
while X remains a viable candidate. 

The annihilation rate of particles of type j with particles of type i is r«"™'*( j") — 
Tii ni{T) < afp'^^'^Vij >, where < ... > indicates a thermal average and Vij is the 
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relative velocity. As the Universe cools, the densities of all the particle species i 
decrease and eventually the rate of even the most important annihilation reaction 
falls below the expansion rate of the Universe. The temperature at which this 
occurs is called the freezeout temperature 7} and can be calculated by solving 
Boltzmann equations. For the freezeout of X annihilation the Boltzmann equation 
can be written as 

X dYx _ T,,rC <av> f \ 

dx H{T) v^l^ 7 ' 

where x — m^/T, Yy = nx/s and s is the entropy of the Universe. Notice 
that dVx/dx goes to zero (or Yx stays constant, corresponding to freezeout) 
when r^™*^(Tjf) <^ H{T). Therefore the freezeout temperature can be estimated 
roughly from a condition Tf-'^'^{Tj) = H{Tj) = 1.66^ Tf/Mpi, where is 
the effective number of relativistic degrees of freedom [28]. Between a few MeV 
and the QCD phase transition only neutrinos, and 7's are in equilibrium and 
= 10.75. Above the QCD phase transition which is within the range 100 to 
200 MeV, light quarks and antiquarks {q, q) and /i"^ are also relativistic species 
in equilibrium, giving g^ — 56.25. The equilibrium density at freeze out temper- 
ature, nj{Tj), is a good estimate of the rehc abundance of the jih species [28]. 
A key element of baryon-number sequestration is that self-annihilation cannot be 
important for maintaining equilibrium prior to freeze out. This is easily satisfied 
as long as cr^"^ is not much greater than cr^" , since at freezeout in "X4" scenario. 

Given mx, Bx and gx (the number of degrees of freedom of the X particle) 
and associated densities n^x^x}-, the temperature Tx at which X's must freeze out 
of thermal equilibrium satisfies: 

r'-x nx nx ^^^Vv^^^ ^ 10.75 r/nilQ-^" g) 



Tlx 2C(3)(27r)3/2 3.91 By 
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where Xx = tnx/Tx- is the factor by which ^ increases above annihilation. 
The equation for X freezeout is the same, with (1 — e) (1 — e)/e. Freezeout 
parameters for our specific models, the H dibaryon and the X, are given in Table 
2.1; a = {a^^^v)/{v) is averaged over the relevant distribution of cm. kinetic 
energies, thermal at ~ 100 MeV for freezeout. 

If Xs interact weakly with nucleons, standard WIMP searches constrain the 
low energy scattering cross section aoM = (c^jv ^'^xn)/{^ + However if the 
X is a hadron, multiple scattering in the earth or atmosphere above the detector 
can cause a significant fraction to reflect or be degraded to below threshold energy 
before reaching a deep underground detector. Scattering also greatly enhances 
DM capture by Earth, since only a small fraction of the halo velocities are less 
than v^g^ — 11 km/s. Table I gives the total fluxes and the factor /cap by which 
the flux of captured X's is lower, for the two scenarios. The capturing rate is 
obtained using code by Edsjo et al. [29] which calculates the velocity distribution 
of weakly interacting dark matter at the Earth taking into account gravitational 
diffusion by Sun, Jupiter and Venus. For the H dibaryon these are the result 
of integrating the conservative halo velocity distribution [30]. A comprehensive 
reanalysis of DM cross section limits including the effect of multiple scattering 
is given in thesis section 4.1 and ref. [30]. A window in the DM exclusion was 
discovered for mx ^2.4 GeV and a^M ~ 0.3 — 1 jjh; otherwise, if the DM mass 
<5 GeV, aoM must be < lO-^^cm^, [30]. 

Since cr^x,x}N is negligible compared to a^N and the X, X do not bind to 
nuclei [31], nucleosynthesis works the same in these scenarios as with standard 
CDM. Primordial light element abundances constrain the nucleon - not baryon - 
to photon ratio! 

The CPT theorem requires that af^'^ + a"^""" = ^ + (7|'°-^''°. Therefore 
a non-trivial consistency condition in this scenario is 
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Table 2.1: Required freezeout temperatures and annihilation cross sections at 
freezeout, and captured DM flux in Earth, in two models; (T_42 = (T/(10~^^cm^). 



Model 


Tx MeV 


Tx MeV 


a^"^ cm^ 


cr^° cm^ 


R 

-"-cap ^5 


H, H 


86.3 


84.5 


2.2 10-^1 


2.8 10-^1 


3.8 X 10^3 


^4 


180 


159 


3.3 10-^^ 


3.7 10-^5 


1.6 X 10^V_42 



CTx -CTjc < 

The value of the LHS needed for B-sequestration from Table I is compatible with 
the upper limits on the RHS from DM searches, and q-^^"^'^^ > cr^, so no flne- 
tuning is required to satisfy CPT. 

Further in the text we focus on specific DM candidates with baryon number 
and experimental constraints that can be placed on such particles. 
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Chapter 3 

H dibaryon: Dark Matter 
candidate within QCD? 

This Chapter is organized as follows: in §3.1 we review properties of the di-baryon; 
in §3.2 we focus on the H which is tightly bound and therefore is light and compact 
object; we review the experiments relevant for the existence of tightly bound H 
in §3.3; we calculate nuclear transitions of the H in §3.4 and binding of the H 
to nuclei in §3.5 and set bounds on parameters from the experiments; in §3.6 
we Summarize the properties of the H which would be consistent with current 
existence experiments. 

3.1 H history and properties 

The H particle is a (udsuds) flavor singlet dibaryon with strangeness -S" = —2, 
charge Q = and spin-isospin-parity = 0"*". In 1977 Jaffe calculated its mass 
[32] to be about 2150 MeV in the MIT-bag model and thus predicted it would be 
a strong-interaction-stable bound state, since decay to two A particles would not 
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be kinematically allowed. The basic mechanism which is expected to give large 
attractive force between quarks in the H is the color magnetic interaction. The 
contribution to the mass from lowest-order gluon exchange is proportional to 



where ai is spin and color vector of the ith quark, and M{miR, rrijR) mea- 
sures the interaction strength. For color singlet hadrons containing quarks and no 
antiquarks [32], 



where N is total number of quarks, J is their angular momentum and Ce is Casimir 
operator for the color-spin representation of quarks, SU{6). We can see that the 
lightest dibaryons will be those in which the quarks are in the color-spin repre- 
sentation with the largest value of Casimir operator. Large values of the Casimir 
operator are associated with symmetric color-spin representation. Antisymmetry 
requires flavor representation to be asymmetric. Calculation shows that only flavor 
singlet representation is light enough to be a stable state. 

This result raised high theoretical and experimental attention. The mass was 
later computed using different models such as Skyrme, quark cluster models, lat- 
tice QCD with values ranging between 1.5 and 2.2 GeV. Such a wide range of 
predictions makes clear contrast to the success in reproducing the mass spectrum 
of mesons and baryons using the above methods. 

The H searches have been done using different production mechanisms, we will 
comment here on the most common approaches 

• H production via {K~,K'^) reaction: In BNL E885 experiment, [33], C12 
is used as a target. The subsequent reactions: K~ + p —>■ + S~ and 
(H~74)atom H -\- X was expected to produce the H] 




(3.1) 




(3.2) 
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• Production through Heavy Ion collision: In BNL E896 experiment, [34], a 
Au + Au collision was used to produce the H that could be identified by the 
anticipated decay H — > TTp — > mr~p; 

• p- nucleus annihilation reaction: In the experiment by Condo et al., [35], 
they used the reaction p + A^H + X to produce the H, and looked for its 
decay through if — > E~ + p channel. 

Other experiments can be found in ref. [36]. 

Experiments were guided by theoretical predictions for H production and decay 
lifetimes. The most remarkable contribution to theoretical predictions came from 
the works of Aerts and Dover [37-39] whose calculations have been the basis for 
understanding of the Brookhaven experiments BNL E813 and E836 with respect 
to the formation of the H. However, theoretical predictions may depend critically 
on the wave function of the H dibaryon. 

Experiments so far did not confirm the existence of the H particle but they 
put bounds on the production cross section for particular mass values, see [36] for 
a detailed review. An underlying assumption has generally been that the H is not 
deeply bound. In our work we are particularly interested in the possibility that 
the H is tightly bound and that it has a mass less than thn + ttia. In that case, 
as we shall see, its lifetime can be longer than the age of the Universe. 

3.2 Tightly bound H 

Several lines of reasoning suggest the possibility that the H could be a tightly 
bound state, with mass lower than mN + tua — 2053 MeV. We digress here briefiy 
to review this motivation, put forth in [40] . The first line of motivation starts from 
the observation that the properties of the | baryon resonance A(1405) and its 
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spin I partner A(1520) are nicely explained if these are assumed to be "hybrid 
baryons": bound states of a gluon with three quarks in a color -octet, flavor- 
singlet state denoted {uds)8- If we adopt the hybrid baryon interpretation of the 
A(1405) and A(1520), the similarity of their masses and glueball mass (~ 1.5 GeV) 
suggests that the color singlet bound state of two {uds)s^s, which would be the H, 
might also have a mass close to 1.5 GcV. A second line of reasoning starts from 
the observation that instantons produce a strong attraction in the scalar diquark 
channel, explaining the observed diquark structure of the nucleon. The H has a 
color-flavor-spin structure which permits a large proportion of the quark pairwise 
interactions to be in this highly attractive spin -0, color 3 channel, again suggesting 
the possibility of tightly bound H. Indeed, ref. [41] reported an instanton-gas 
estimate giving niH — 1780 MeV. If the H is tightly bound, it is expected to be 
spatially compact. Hadron sizes vary considerably, for a number of reasons. The 
nucleon is signiflcantly larger than the pion, with charge radius tat = 0.87 fm 
compared to = 0.67 fm [5]. Lattice and instanton- liquid studies qualitatively 
account for this diversity and further predict that the scalar glueball is even more 
tightly bound: ~ 0.2 fm [42,43]. If the analogy suggested in ref. [44] between 
H, Ai405 and glueball is correct, it would suggest rn ~ rG<l/4 vn- The above 
size relationships make sense: the nucleon's large size is due to the low mass of the 
pion which forms an extended cloud around it, while the H and glueball do not 
couple to pions, due to parity and flavor conservation, are thus are small compared 
to the nucleon. 

Lattice QCD efforts to determine the H mass have been unable to obtain clear 
evidence for a bound H. The discussion above suggests that inadequate spatial 

resolution may be a problem due to the small size of the H. Several lattice calcula- 
tions [45-48] use sufficiently small spacing but they use quenching approximation, 
which does not properly reproduce instanton effects. This is a serious deficiency 
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since the instanton liquid results of ref. [41] indicate that instanton contributions 
are crucial for the binding in the scalar diquark channel. In the absence of an 
unquenched, high-resolution lattice QCD calculation capable of a reliable deter- 
mination of the H mass and size, we will consider all values of uih and take 
'"h/^^tv = 1// as a parameter, with / in the range 2-6. 

Based on the fact that the H interacts with ordinary matter only if it has non- 
vanishing color charge radius (since it is spin and isospin zero neutral particle) 
ref. [40] estimates cross section of the H with ordinary matter to be of the order 
10~^'^ mb. 

Based on the assumption of light and tightly bound H motivated by Farrar in 
[40] we examine the viability of this model using current experimental constraints. 
My work has focused on the study of several processes involving the H which are 
phenomenologically important if it exists: 

• whether it binds to nuclei, in which case exotic isotope searches place strin- 
gent constraints; 

• nucleon transitions of the H - conversion of two A's in a doubly-strange 
hypernucleus to an H which is directly tested in experiments, decay of the 
H to two baryons, and — if the H is light enough — conversion of two nucleons 
in a nucleus to an if. 

The experimental constraints important for the existence of the light H and the 
if as a DM candidate are outlined below. For more experimental constraints on 
tightly bound H, see [40]. 
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3.3 The Existence of the H - Experimental con- 
straints 



3.3.1 Double A hyper-nucleus detection 

One of the ways to produce and study the H is through the production of double 
hypernuclei. A double A hypernucleus formed in an experiment usually through 
the {K~^K^) interaction with the target, is expected to decay strongly into the 
H and the original nucleus. If single A decay from a double A hypernucleus is 
observed, it means that either the mass of the H should be heavier than the mass 
of the two A's minus the binding energy, or that the decay of two A's to an H 
proceeds on a longer time scale than the A weak decay. There are five experiments 
which have reported positive results in the search for single A decays from double 
A hypernuclei. The three early emulsion based experiments [49-51] suffer from 
ambiguities in the particle identification, and therefore we do not consider them 
further. In the latest emulsion experiment at KEK [52] , an event has been observed 
which is interpreted with good confidence as the sequential decay of He^y^ emitted 
from a S~ hyperon nuclear capture at rest. The binding energy of the double A 
systems is obtained in this experiment to be Saa = 1.01 ± 0.2 MeV, in significant 
disagreement with the results of previous emulsion experiments, finding Saa ~ 4.5 
MeV. 

The third experiment at BNL [53] was not an emulsion experiment. After 
the {K~^K^) reaction on a Bc^ target produced S=-2 nuclei it detected tt pairs 
coming from the same vertex at the target. Each pion in a pair indicates one unit 
of strangeness change from the (presumably) di-A system. Observed peaks in the 
two pion spectrum have been interpreted as corresponding to two kinds of decay 
events. The pion kinetic energies in those peaks are (114,133) MeV and (104,114) 
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MeV. The first peak can be understood as two independent single A decays from 
A A nuclei. The energies of the second peak do not correspond to known single 
A decay energies in hyper-nuclei of interest. The proposed explanation [53] is 
that they are pions from the decay of the double A system, through a specific 
He resonance. The required resonance has not yet been observed experimentally, 
but its existence is considered plausible. This experiment does not suffer from low 
statistics or inherent ambiguities, and one of the measured peaks in the two pion 
spectrum suggests observation of consecutive weak decays of a double A hyper- 
nucleus. The binding energy of the double A system could not be determined 
in this experiment. 

The KEK and BNL experiments are generally accepted to demonstrate quite 
conclusively, in two different techniques, the observation of A decays from double 
A hypernuclei. Therefore the formation of the in a double A hypernucleus does 
not proceed, at least not on a time scale faster than the A lifetime, i.e., ta/^j^-^a'^x 
cannot be much less than Ri 10~-^°s. (To give a more precise limit on ta^a^a^x 
requires a detailed analysis by the experimental teams, taking into account the 
number of hypernuclei produced, the number of observed A decays, the acceptance, 
and so on.) This experiment is considered leading evidence against the existence 
of the H di-baryon. As will be seen below, this constraint is readily satisfied if the 
H is compact: r// ^ 1/2 r^r or less, depending on the nuclear wave function. 

3.3.2 Stability of nuclei 

This subsection derives constraints for a stable iJ, where itlh < 2mjv- In that case 
the H would be an absolutely stable form of matter and nuclei would generally 
be unstable toward decay to the H. There are a number of possible reactions 
by which two nucleons can convert to an in a nucleus if that is kinematically 
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allowed {mn ^'2mN^). The initial nucleons are most likely to be pn or nn in 
a relative s-wave, because in other cases the Coulomb barrier or relative orbital 
angular momentum suppresses the overlap of the nucleons at short distances which 
is necessary to produce the H. If mn^'^rnN — rriT^ — 11 AQ MeV, the final state 
can be Ht:'^ or Htt^. If tt production is not allowed, for mn ^ 1740 MeV, the 
most important reactions are pn — > He^Ve or the radiative-doubly-weak reaction 
nn — > H^. 

The best experiments to place a limit on the stability of nuclei are proton 
decay experiments. Super Kamiokande (SuperK), can place the most stringent 
constraint due to its large mass. It is a water Cerenkov detector with a 22.5 
kiloton fiducial mass, corresponding to 8 10^^ oxygen nuclei. SuperK is sensitive 
to proton decay events in over 40 proton decay channels [54]. Since the signatures 
for the transition of two nucleons to the H are substantially different from the 
monitored transitions, a specific analysis by SuperK is needed to place a limit. We 
will discuss the order-of-magnitude of the limits which can be anticipated. 

Detection is easiest if the H is light enough to be produced with a tt"^ or . 
The efficiency of SuperK to detect neutral pions, in the energy range of interest 
(KE ~ — 300 MeV), is around 70 percent. In the case that a tt"^ is emitted, 
it can charge exchange to a 7r° within the detector, or be directly detected as 
a non-showering muon-like particle with similar efficiency. More difficult is the 
most interesting mass range m//> 1740 MeV, for which the dominant channel 
pn — > He^v gives an electron with E ~ (2m jv — mij)/2< 70 MeV. The channel 
nn — >H7, whose rate is smaller by a factor of order a, would give a monochromatic 
photon with energy (2m7v — rnn) ^ 100 MeV. 

We can estimate SupcrK's probable sensitivity as follows. The ultimate back- 

"'^ Throughout, wc use this shorthand for the more precise incquahty niH < tua — tua' — ttix 
where nix is the minimum invariant mass of the final decay products. 
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ground comes primarily from atmospheric neutrino interactions, 

uN ^ N'{e,^i), i^N ^ N'{e,fi) +n7: anduN ^ uN' + n7r, (3.3) 

for which the event rate is about 100 per kton-yr. Without a strikingly distinct 
signature, it would be difficult to detect a signal rate significantly smaller than 
this, which would imply SuperK might be able to achieve a sensitivity of order 
TAj^ff-^A'^x ^ fewlO^^ yr. Since the H production signature is not more favorable 
than the signatures for proton decay, the SuperK limit on taj^,^-*a'jjX can at best be 
O.lTp, where 0.1 is the ratio of Oxygen nuclei to protons in water. Detailed study 
of the spectrum of the background is needed to make a more precise statement. 
We can get a lower limit on the SuperK lifetime limit by noting that the SuperK 
trigger rate is a few Hz [54], putting an immediate limit To_>ij+x ^ few 10^^ yr, 
assuming the decays trigger SuperK. 

SuperK limits will apply to specific decay channels, but other experiments 
potentially establish limits on the rate at which nucleons in a nucleus convert to 
an H which are independent of the H production reaction. These experiments 
place weaker constraints on this rate due to their smaller size, but they are of 
interest because in principle they measure the stability of nuclei directly. Among 
those cited in ref. [55], only the experiment by Flerov et. al. [56] could in principle 
be sensitive to transitions of two nucleons to the H. It searched for decay products 
from Th^^^, above the Th natural decay mode background of 4.7 McV a particles, 
emitted at the rate r„ = 0.7 10"^^^"^ Cuts to remove the severe background of 
4.7 MeV a's may or may not remove events with production of an H. Unfortunately 
ref. [56] does not discuss these cuts or the experimental sensitivity in detail. An 
attempt to correspond with the experimental group, to determine whether their 
results are applicable to the H, was unsuccessful. If applicable, it would establish 
that the lifetime TTh232_>^^;5f > lO^-*^ yr. 
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Better channel-independent limits on N and A^A^ decays in nuclei have been 
established recently, as summarized in ref. [57]. Among them, searches for the 
radioactive decay of isotopes created as a result of NN decays of a parent nucleus 
yield the most stringent constraints. This method was first exploited in the DAMA 
hquid Xe detector [58]. BOREXINO has recently improved these results [57] using 
their prototype detector, the Counting Test Facility (CTF) with parent nuclei 
C^^, C^^ and O^^. The signal in these experiments is the beta and gamma radiation 
in a specified energy range associated with deexcitation of a daughter nucleus 
created by decay of outer-shell nucleons in the parent nucleus. They obtain the 
hmits Tpp > 5 10^^ yr and t„„ > 4.9 10^^ yr. However H production requires overlap 
of the nucleon wavefunctions at short distances and is therefore suppressed for 
outer shell nucleons, severely reducing the utility of these limits. Since the Super K 
limits will probably be much better, we do not attempt to estimate the degree of 
suppression at this time. 

Another approach could be useful if for some reason the direct SuperK search is 
foiled. Ref. [59] places a limit on the lifetime of a bound neutron, r„ > 4.9 10^^ yr, 
by searching for 7's with energy = 19 — 50 MeV in the Kamiokande detector. 
The idea is that after the decay of a neutron in oxygen the de-excitation of O^^ 
proceeds by emission of 7's in the given energy range. The background is especially 
low for 7's of these energies, since atmospheric neutrino events produce 7's above 
100 MeV. In our case, some of the photons in the de-excitation process after 
conversion of to H, would be expected to fall in this energy window. 

In §3.4 we calculate nuclear transition rates of a tightly bound H in order to 
find constraints on the H size and mass, primarily focusing on the constraints set 
by SuperK experiment. 
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3.3.3 Experimental constraints on the H binding 

If the H binds to nuclei and if it is present with DM abundance, it should be 
present in nuclei on Earth and experiments searching for anomalous mass isotopes 
would be sensitive to its existence. Accelerator mass spectroscopy (AMS) experi- 
ments generally have high sensitivity to anomalous isotopes, limiting the fraction 
of anomalous isotopes to 10"^^ depending on the element. We discuss binding of 
the H to heavy and to hght isotopes separately. 

The H will bind more readily to heavy nuclei than to light ones because their 
potential well is wider. However, searches for exotic particles bound to heavy nuclei 
are limited to the search for charged particles in Fe [60] and to the experiment 
by Javorsek et al. [61] on Fe and Au. The experiment by Javorsek searched for 
anomalous Au and Fe nuclei with Mx in the range 200 to 350 atomic mass units 
u. Since the mass of Au is 197 u, this experiment is sensitive to the detection of an 
exotic particle with mass Mx > 3 u= 2.94 GeV and is not sensitive to the tightly 
bound H. 

A summary of limits from various experiments on the concentrations of exotic 
isotopes of light nuclei is given in [62] . Only the measurements on hydrogen [63] and 
helium [64] nuclei are of interest here because they are sensitive to the presence 
of a light exotic particle with a mass of Mx ~ 1 GeV. It is very improbable 
that the H binds to hydrogen, since the A does not bind to hydrogen in spite of 
having attractive contributions to the potential not shared by the H, e.g., from 
r] and rj' exchange. Thus we consider only the limit from helium. The limit on 
the concentration ratio of exotic to non-exotic isotopes for helium comes from 
the measurements of Klein, Middleton and Stevens who quote an upper limit of 
^ < 2 X 10"^^ and ^ < 2 X 10"^^ for primordial He [65]. Whether these 
constraints rule out the H depends on the H coupling to nucleus. 
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In §3.5 we calculate the binding of the H, or more generally any flavor singlet, 
to nuclei and find the values of coupling which are allowed from the existence of 
the H. As we will see, the allowed couplings coincide with the values expected 
from the particle physics arguments. 

3.4 Nucleon and nuclear transitions 

of the H dibaryon - Rates estimate 

In this section we study several processes involving the H which are phenomeno- 
logically important if it exists, [66]: conversion of two A's in a doubly-strange 
hypernucleus to an H (§3.3.1), decay of the H to two baryons, and — if the H 
is light enough — conversion of two nucleons in a nucleus to an H (§3.3.2). The 
amplitudes for these processes depend on the spatial wavefunction overlap of two 
baryons and an H. We are particularly interested in the possibility that the H is 
tightly bound and that it has a mass less than ttin + rn\ because then, as we shall 
see, the H is long-hved, with a lifetime which can be longer than the age of the 
Universe. 

To estimate the rates for these processes requires calculating the overlap of 
initial and final quark wavef unctions. We model that overlap using an Isgur-Karl 
harmonic oscillator model for the baryons and H, and the Bethe-Goldstone and 
Miller-Spencer wavefunctions for the nucleus. The results depend on tn/th and 
the nuclear hard core radius. 

We also calculate the lifetime of the H taking similar approach to the overlap 
calculation. We do it in three qualitatively distinct mass ranges, under the as- 
sumption that the conditions to satisfy the constraints from double-A hypernuclei 
are met. The ranges are 
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• tuh < rriN + tua, in which H decay is a doubly- weak A^* = 2 process, 

• niN + TTiA < rriH < 2mA, in which the H can decay by a normal weak 
interaction, and 

• ruH > '^TTiA, in which the H is strong-interaction unstable. 

The H hfetime in these ranges is greater than or of order 10^ years, ~ 10 sec, and 
~ 10~^^ sec, respectively. 

Finally, if mjy <2mAr, nuclei are unstable and AS* = —2 weak decays convert 
two nucleons to an H. In this case the stability of nuclei is a more stringent 
constraint than the double-A hypernuclear observations, but our results of the 
next subsection show that nuclear stability bounds can also be satisfied if the H 
is sufficiently compact: ^ 1/4 rN depending on mass and nuclear hard core 
radius. This option is vulnerable to experimental exclusion by SuperK. 

In order to calculate transition rates we factor transition amplitudes into an 
amplitude describing an i7-baryon-baryon wavefunction overlap times a weak inter- 
action transition amplitude between strange and non-strange baryons. In subsec- 
tion 3.4.1 we setup the theoretical apparatus to calculate the wavefunction overlap 
between H and two baryons. We determine the weak interaction matrix elements 
phenomenologically in subsection 3.4.2. Nuclear decay rates are computed in sub- 
section 3.4.3 while lifetime of the H for various values of ttih is found in 3.4.4. The 
results are reviewed and conclusions are summarized in section 3.6. 

3.4.1 Overlap of H and two baryons 

We wish to calculate the amplitudes for a variety of processes, some of which 
require one or more weak interactions to change strange quarks into light quarks. 
By working in pole approximation, we factor the problem into an if-baryon-baryon 
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wavefunction overlap times a weak interaction matrix element between strange and 
non-strange baryons, which will be estimated in the next section. For instance, 
the matrix element for the transition of two nucleons in a nucleus A to an H and 
nucleus A', Ajvat -^h-^i calculated in the AA pole approximation, as the 
product of matrix elements for two subprocesses: a transition matrix element for 
formation of the H from the A A system in the nucleus, |A1|{AA}^/f x, times the 
amplitude for a weak doubly-strangeness-changing transition, IA^Iatat-^aa: 



We ignore mass differences between light and strange quarks and thus the spatial 
wavefunctions of all octet baryons are the same. In this section we are concerned 
with the dynamics of the process and we suppress spin-flavor indices. 

Isgur-Karl Model and generalization to the H 

The Isgur-Karl (IK) non-relativistic harmonic oscillator quark model [67-69] was 
designed to reproduce the masses of the observed resonances and it has proved to 
be successful in calculating baryon decay rates [68]. In the IK model, the quarks 
in a baryon are described by the Hamiltonian 



where we have neglected constituent quark mass differences. The wave function of 
baryons can then be written in terms of the relative positions of quarks and the 
center of mass motion is factored out. The relative wave function in this model 



\M\a^A'^X — \M\{KK}^H X IA^IatAT^AA- 



(3.4) 



(3.5) 



is [68, 69] 



*b(^i, ^2, nj) = Nb exp 




(3.6) 
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where Nb is the normahzation factor, as = , ^ ^ = Vskm, and < Tr > is the 
baryon mean charge radius squared. Changing variables to 

^ fi - r1 - ri + f2- 2r3 

""^•^^ — 

reduces the wave function to two independent harmonic oscillators. In the ground 
state 



(3.8) 



One of the deficiencies of the IK model is that the value of the parame- 
ter needed to reproduce the mass sphttings of lowest lying and I"*" baryons, 
as = 0.406 GeV, corresponds to a mean charge radius squared for the proton of 
^y < rj^ > = ^ = 0.49 fm. This is distinctly smaller than the experimental value 
of 0.87 fm. Our results depend on the choice of ckb and therefore we also report 
results using as — 0.221 GeV which reproduces the observed charge radius at the 
expense of the mass-splittings. 

Another concern is the applicability of the non-relativistic IK model in de- 
scribing quark systems, especially in the case of the tightly bound H. With 
fn/fN — 1//, the quark momenta in the H are ~ / times higher than in the 
nucleon, which makes the non-relativistic approach more questionable than in the 
case of nucleons. Nevertheless we adopt the IK model because it offers a tractable 
way of obtaining a quantitative estimate of the effect of the small size of the H 
on the transition rate, and there is no other alternative available at this time. For 
comparison, it would be very interesting to have a Skyrme model calculation of 
the overlap of an H with two baryons. 

We fix the wave function for the H particle starting from the same Hamiltonian 
(3.5), but generalized to a six quark system. For the relative motion part this gives 



"^H — Nff exp 



i<j 



(3.9) 
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The space part of the matrix element of interest, {A'jj\A\\), is given by the integral 

6 

/ n d'rinih 2, 3)^1(4, 5, 6)^^(1, 2, 3, 4, 5, 6). (3.10) 

i=i 

Therefore it is useful to choose variables for the H wavefunction as follows, replac- 
ing 

n, r2, fa, n, fg, fe ^, A", p*, A^ a, Rcm, (3.11) 

where p***^^^ and A"*^^^ are defined as in eq. (3.7), with a{b) referring to coordinates 
1, 2, 3 (4, 5, 6). (Since we are ignoring the flavor-spin part of the wavefunction, we 
can consider the six quarks as distinguishable and not worry about fermi statistics 
at this stage.) We also define the center-of-mass position and the separation, a, 
between initial baryons a and b: 

pa _|_ r>b 

3 _ ^CM "T -"•CM pa p6 fi ^0\ 

Using these variables, the H ground state wave function becomes 

o\ 3/4 / N 15/2 
3 \ I C^H 



'n ^ 

r a% , ^2 -2 -2 ^2 3 

X exp[ — ^(p« +A« +A'' +2«)]- 



As for the 3-quark system, an = — ^ 



Nucleeir Wavefunction 



We will use two different wavefunctions to describe two A's or nucleons in a nu- 
cleus, in order to study the model dependence of our results and to elucidate the 
importance of different aspects of the nuclear wavefunction. A commonly used 
wavefunction is the Miller- Spencer (MS) wavefunction [70]: 

V'M5 = l-exp-^i«'(l-C2a2), (3.14) 
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with the canonical parameter choices Ci = 1.1 fm"^ and C2 = 0.68 fm~^. It must 
be emphasized that at the short distances relevant for this calculation, the form 
and magnitude of the MS wavefunction are not constrained experimentally and 
rather are chosen to give a good fit to long-distance physics with a simple func- 
tional form. The other wavefunction we use (BBG) is a solution of the Bruecker- 
Bethe-Goldstone equation describing the interaction of a pair of fermions in an 
independent pair approximation; see, e.g., [71]. It is useful because we can ex- 
plicitly explore the sensitivity of the result to the unknown short-distance nuclear 
physics by varying the hard-core radius. 

The BBG wave function is obtained as follows. The solution of the Schrodinger 
equation for two fermions in the Fermi sea interacting through a potential v{xi, X2) 
takes the form 

^(1,2)^^6^'''^'^^ ^l^ia), (3.15) 
V v 

where Rcm and a are defined as in (3.12). The first factor contains the center-of- 
mass motion and the second is the internal wave function of the interacting pair. 
ip{a) is a solution of the Bethe-Goldstone equation (eq. (36.15) in [71]) which is 
simply the Schrodinger equation for two interacting fermions in a Fermi gas, where 
the Pauli principle forbids the appearance of intermediate states that are already 
occupied by other fermions. Both wave functions are normalized so that the space 
integral of the modulus-squared of the wave function equals one. In the application 
of this equation to nuclear matter, the interaction of each particle from the pair 
with all particles in the nucleus through an effective single particle potential is 
included, in the independent pair approximation known as Bruecker theory (see 
eq. (41.1) and (41.5) in [71]). 

We are interested in s-wave solutions to the Bethe-Goldstone equation since 
they are the ones that penetrate to small relative distances. Following [71], an 
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s-wave solution of the internal wave function is sought in the form 

i^ia) - (3.16) 

(X 

which simphfies the Bethe-Goldstone equation to 

— + k'^j u{a) = v{a)u{a) - J x{a, y)v{y)u{y)dy (3.17) 

where v{a) is the single particle potential in the effective- mass approximation, and 
the kernel x(a, y) is given by 



x(a,y) = - 

TT 



smkp{a — y) smkp{a + y) 



(3.18) 



a — y a, + y 

where kp is the Fermi wavenumber. For the interaction potential between two nu- 
cleons in a nucleus we choose a hard core potential for the following reasons. The 
two particle potential in a nucleus is poorly known at short distances. Measure- 
ments (the observed deuteron form factors, the sums of longitudinal response of 
light nuclei,...) only constrain two-nucleon potentials and the wave functions they 
predict at internucleon distances larger than 0.7 fm [72]. The Bethe-Goldstone 
equation can be solved analytically when a hard-core potential is used. While the 
hard-core form is surely only approximate, it is useful for our purposes because it 
enables us to isolate the sensitivity of the results to the short-distance behavior of 
the wavefunction. We stress again, that more "realistic" wavef unctions, including 
the MS wave function, are in fact not experimentally constrained for distances be- 
low 0.7 fm. Rather, their form at short distance is chosen for technical convenience 
or aesthetics. 

Using the hard core potential, the s-wave EG wavefunction is 

Nbg^ for a>^ 
*BG(a) = , (3.19) 

for a < 7^ 
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with 



Nbg = , (3.20) 



u{a) 



47r da 



where ^ is the hard core radius and R{A) — 1.07^4^/^ is the radius of a nucleus 
with mass number A. Expressions for u can be found in [71], eq. (41.31). The 
normahzation factor Nbg is fixed setting the integral of IV'bgP over the volume 
of the nucleus equal to one. The function u vanishes at the hard core surface 
by construction and then rapidly approaches the unperturbed value, crossing over 
that value at the so called "healing distance" . At large relative distances and when 
the size of the normalization volume is large compared to the hard core radius, 
u{a)/a approaches a plane wave and the normalization factor Nbg (3.20) reduces 
to the value 1 / VHox, as 

i^BGia) ^Nbg^ ^ e'^\ (3.21) 

Overlap Calculation 

The non-relativistic transition matrix element for a transition AA — > H inside a 
nucleus is given by (suppressing spin and flavor) 

T{KK}^H = 2m5{E) I d^a d^RcM \{ d^p'd^X' 

i=a,b 

X ^H^l < i^nuc e^(^«-^^^)^^^, (3.22) 



where 5(E) = 6iEH - ^aa), i^f = i^fif^'^'M), and ^nuc = ^nuc{a) is the 
relative wavefunction function of the two A's in the nucleus. The notation {AA} is 
a reminder that the A's are in a nucleus. The plane waves of the external particles 
contain normalization factors 1 / \/V and these volume elements cancel with volume 
factors associated with the final and initial phase space when calculating decay 
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rates. The integration over the center of mass position of the system gives a 3- 
dimensional momentum delta function and we can rewrite the transition matrix 
element as 

T{!,K}^H = {2T:)\5\kf - ki) M{AA}^H, (3.23) 

where \ A4\{aa}-,h is the integral over the remaining internal coordinates in eq. (3.22) 
In the case of pion or lepton emission, plane waves of the emitted particles should 
be included in the integrand. For brevity we use here the zero momentum transfer, 

— * 

k = approximation, which we have checked holds with good accuracy; this is not 
surprising since typical momenta are ^0.3 GeV. 

Inserting the IK and BBG wavefunctions and performing the Gaussian integrals 
analytically, the overlap of the space wave functions becomes 

X Nbg / d^a^^e-i'^'i^ 

where the factor l/-\/4 comes from the probability that two nucleons are in a 
relative s-wave, and / is the previously-introduced ratio of nucleon to H radius; 
otH = f CUB- Since Nbg has dimensions the spatial overlap A4{aa}-^h is 

a dimensionless quantity, characterized by the ratio /, the Isgur-Karl oscillator 
parameter a^, and the value of the hard core radius. Fig. 3.1 shows |A^|{AA}-»if 
calculated for oxygen nuclei, versus the hard-core radius, for a range of values of 
/, using the standard value of q;b = 0.406 GeV for the IK model [69] and also 
aB = 0.221 GeV for comparison. 

Fig. 3.1 shows that, with the BBG wavefunction, the overlap is severely sup- 
pressed and that the degree of suppression is very sensitive to the core radius. 
This confirms that the physics we are investigating depends on the behavior of the 
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Figure 3.1: Logio of \M.W/^^h versus hard core radius in fm, for ratio / = Rn/Rh 
and two values of the Isgur-Karl oscillator parameter: as = 0.406 GeV (thick 
lines) and = 0.221 GeV (thin lines). 

nuclear wavefunction at distances at which it is not directly constrained experi- 
mentally. Fig. 3.2 shows a comparison of the overlap using the Miller Spencer and 
BBG nuclear wavef unctions, as a function of the size of the H. One sees that the 
spatial overlap is strongly suppressed with both wavef unctions, although quanti- 
tatively the degree of suppression differs. We cannot readily study the sensitivity 
to the functional form of the baryonic wavefunctions, as there is no well-motivated 
analytic form we could use to do this calculation other than the IK wavefunction. 
However by comparing the extreme choices of parameter cts in the IK wavefunc- 
tion, also shown in Figs. 3.1 and 3.2, we explore the sensitivity of the spatial 
overlap to the shape of the hadronic wavefunctions. Fortunately, we will be able 
to use additional experimental information to constrain the wavefunction overlap 
so that our key predictions are insensitive to the overlap uncertainty. 
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Figure 3.2: Logio of |-^Iaa-*h versus ratio / = an/oiN, calculated with BBG wave 
function with core radius 0.4 and 0.5 fm, and with the MS wave function. Thick 
(thin) lines are for = 0.406 GeV {ptB = 0.221 GeV) in the IK wavefunction. 

3.4.2 Weak Interaction Matrix Elements 

Transition of a two nucleon system to off-shell AA requires two strangeness chang- 
ing weak reactions. Possible /S.S = 1 sub-processes to consider are a weak transition 
with emission of a pion or lepton pair and an internal weak transition. These are 
illustrated in Fig. 3.3 for a three quark system. We estimate the amplitude for 
each of the sub-processes and calculate the overall matrix element for transition 
to the AA system as a product of the sub-process amplitudes. 



The matrix element for weak pion emission is estimated from the A — > iVvr 
rate: 



l-^li^iv.= 7^ ^-^ = 0.8x10-- GeV^ (3.25) 

By crossing symmetry this is equal to the desired in the approximation 

of momentum-independence which should be valid for the small momenta in this 
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Figure 3.3: Some relevant weak transitions for NN — > HX 

application. Analogously, for lepton pair emission we have 

= 7^ - 3 X 10--. (3.26) 

The matrix element for internal conversion, {uds) {udd), is proportional to 
the spatial nucleon wave function when two quarks are at the same point: 

M a^n ~< V'A 53 n - r2 V^iv > — -, (3.27) 

rriq 

where rUq is the quark mass introduced in order to make the 4 point vertex am- 
plitude dimensionless [73]. The expectation value of the delta function can be 
calculated in the harmonic oscillator model to be 

< i^A\S\n - f2)\il^N > = (^) = X 10"' GeV^ (3.28) 



The delta function term can also be inferred phenomenologically in the following 
way, as suggested in [73] . The Fermi spin-spin interaction has a contact character 
depending on (7l(T^/mg5(fi — f2), and therefore the delta function matrix element 
can be determined in terms of electromagnetic or strong hyperfine splitting: 

(mso - ms+) - (m„ - nip) = < - ^2) > (3.29) 

Stt 

mA - = "5^^ < S^in - ra) >, (3.30) 



40 



where rriq is the quark mass, taken to be tjin/^- Using the first form to avoid the 
issue of scale dependence of as leads to a value three times larger than predicted 
by the method used in eq. (3.28), namely: 

< i^Al^^ifi - f2)\ipN > = 1.2 X 10-2 GeV^ (3.31) 

We average the expectation values in eq. (3.28) and eq. (3.31) and adopt 

\M\l^N = X 10-^^ . (3.32) 

In this way we have roughly estimated all the matrix elements for the relevant 
sub-processes based on weak-interaction phenomenology. 

3.4.3 Nuclear decay rates 
Lifetime of doubly-strange nuclei 

The decay rate of a doubly-strange nucleus is: 

where $2 is the two body final phase space factor, defined as in [55], and ttiaa 
is the invariant mass of the A's, 2m\. The factor K contains the transition 
element in spin flavor space. It can be estimated by counting the total number of 
flavor-spin states a uuddss system can occupy, and taking to be the fraction 
of those states which have the correct quantum numbers to form the H. That 
gives K"^ ~ 1/1440, and therefore we write K'^ — (1440 Kiua)'^- Combining these 
factors we obtain the estimate for the formation time of an if in a doubly-strange 
hypernucleus 

^ 3(7) /Vino 10'^' s 
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where the phase space factor was evaluated for mn = 1-8(2) GeV. 

Fig. 3.2 shows |A^|{AA}-*ij the range of / and hard-core radius where its value 
is in the neighborhood of the experimental limits, for the standard choice as — 
0.406 GeV and comparison value ub — 0.221 GeV. In order to suppress r{AAA — > 
A'jjX) sufficiently that some A's in a double-A hypernucleus will decay prior to 
formation of an H, we require I^a^j/ ~ 10^^. If the nucleon hard core potential 
is used, this is satisfied even for relatively large H, e.g., th '^tn/'^.S (rjv/2.1) for a 
hard-core radius 0.4 (0.5) fm and can also be satisfied with the MS wave function 
as can be seen in Fig. 3.2. Thus the observation of single A decay products from 
double-A hypernuclei cannot be taken to exclude the existence of an H with mass 
below 2mA unless it can be demonstrated that the wave function overlap is large 
enough. 

Conversion of AS — nucleus to an H 

If the H is actually stable {tjih < 2mp + 2me) two nucleons in a nucleus may 
convert to an H and cause nuclei to disintegrate. A^A^ HX requires two weak 
reactions. If mH < 1740 MeV two pion emission is allowed and the rate for the 
process Ann ^ ■^'h^'^j approximately 

'-^ - i^^l^ ^3 (3.35) 

V (2mA - J 
where the denominator is introduced to correct the dimensions in a way suggested 
by the AA pole approximation. Since other dimensional parameters relevant to 
this process, e.g., nig — miv/3 or Aqcd, are comparable to 2mA — "m-H and we are 
only aiming for an order-of-magnitude estimate, any of them could equally well be 
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used. The lifetime for nuclear disintegration with two pion emission is thus 

40 K1440 , . 

TA^^^A'^nn ~ TTJT2 (3.36) 

taking ttih — 1-5 GeV in the phase space factor. For the process with one pion 
emission and an internal conversion, our rate estimate is 



r.„„-.,, « K'^^ (3.37) 

X {\M\N^AiT\M\N^A\M\AA^Hf, 



leading to a hfetime, for mn — 1-5 GeV, of 

3 /ti440 . s 

TAnn-^A'^tt ~ TJ-7T2 (3.38) 

If niH ^ 1740 MeV, pion emission is kinematically forbidden and the relevant 
final states are e'^i/ or 7; we now calculate these rates. For the transition Ajvat — 
A'fjeu, the rate is 

X {\M\N^Aeu\M\N^A\M\AA^Hf ■ 

In this case, the nuclear lifetime is 

TA,,-.A',e. ^ I 10^ yr, (3.40) 

taking niH — 1-8 GeV. For Ann A'jjj, the rate is approximately 

2 

X 4>2(|At|5,^A|A1|AA-ff)', 

leading to the lifetime estimate 

TA^^^A'^j ~ 10' yr, (3.42) 

K^i|AA->i? 
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for rriH = 1.8 GeV. 

One sees from Fig. 3.1 that a lifetime bound of ^few 10^^ yr is not a very 
stringent constraint on this scenario if uih is large enough that pion final states 
are not allowed. E.g., with K1440 = 1 the rhs of eq. (3.40) is ^few 10^^ yr, for 
standard aB, hard core radius of 0.45 fm, and ^ 1/5 r^r — in the middle of 
the range expected based on the glueball analogy. If ttih is light enough to permit 
pion production, experimental constraints are much more powerful. We therefore 
conclude that mn ^ 1740 MeV is disfavored and is likely to be excluded, depending 
on how strong limits SuperK can give. 

Table 3.1: The final particles and momenta for nucleon-nucleon transitions to H 
in nuclei. For the 3-body final states marked with *, the momentum given is for 

the configuration with H prodTiced at rest. 



mass 
ruH [GeV] 


final state 
A' H + 


final momenta 

p [MeV] 


partial lifetime 

xK'\M\l^^H [yr] 


1.5 


TT 


318 


2 10"^ 


1.5 


TTTT 


170* 


0.03 


1.8 


eu 


48* 


70 


1.8 


7 


96 


2 10^ 



3.4.4 Lifetime of an Unstable H 

If 2miv ^ niH < mjv + in a, the H is not stable but it proves to be very long lived if 
its wavefunction is compact enough to satisfy the constraints from doubly-strange 
hypernuclei discussed in section 3.4.3. The limits on nuclear stability discussed in 
the previous section do not apply here because nuclear disintegration to an H is 
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not kinematically allowed. 



Wavefunction Overlap 

To calculate the decay rate of the H we start from the transition matrix element 
eq. (3.22). In contrast to the calculation of nuclear conversion rates, the outgoing 
nucleons are asymptotically plane waves. Nonetheless, at short distances their 
repulsive interaction suppresses the relative wavefunction at short distances much 
as in a nucleus. It is instructive to compute the transition amplitude using two 
different approximations. First, we treat the nucleons as plane waves so the spatial 
amplitude is: 

Th^aa = 27riS{EAA - Eh) J J] d^p'd^>^'d^a d^RcM 

i=a,b 

X -0i?V'r V'A e'^^^^^^"^''^'^''''- (3-43) 



The integration over Rcm gives the usual 4D 5 function. Using the Isgur-Karl 
wave function and performing the remaining integrations leading to \J^\h^aa, as 
in eq. 3.23), the amplitude is: 

2/ Vf3\"Ua^Y'' 



X 



/ d^a e ' 2 « 



sy^Y 2/ y _3/2 

3nj \l + f 



The amplitude depends on the size of the H through the factor f — rAf/rn- Note 
that the normalization Nbg in the analogous result eq. 3.24) which comes from 
the Bcthe-Goldstone wavefunction of A's in a nucleus has been replaced in this 
calculation by the plane wave normalization factor l/\/V which cancels with the 
volume factors in the phase space when calculating transition rates. 
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Transition rates calculated using eq. (3.44) provide an upper limit on the true 
rates, because the calculation neglects the repulsion of two nucleons at small dis- 
tances. To estimate the effect of the repulsion between nucleons we again use the 
Bethe-Goldstone solution with the hard core potential. It has the desired prop- 
erties of vanishing inside the hard core radius and rapidly approaching the plane 
wave solution away from the hard core. As noted in section 3.4.1, Nbg 
for a — > oo. Therefore, we can write the transition amplitude as in eq. (3.24), with 
the normalization factor l/\/V canceled with the phase-space volume element: 



\M\ 



H-^AA — 



1 + /V \2j VV^. 
X / d^a^e-t^"'. (3.45) 



Table 3.2: \M\ 

ff-*AA GeV for different values of / (rows) and nuclear wave- 
function (columns), using the standard value asi — 0.406 GeV and the comparison 
value aB2 = 0.221 GeV in the IK wavefunction of the quarks. 





BBG, 0.4 fm 


BBG, 0.5 fm 


MS 














"52 


4 


6 10-1^ 


6 10-8 


7 10-18 


4 10-9 


1 10-8 


8 10-^ 


3 


5 10-9 


3 10-5 


3 10-11 


7 10-6 


2 10-6 


9 10-5 


2 


1 IQ-i 


0.02 


1 10-^ 


0.01 


9 10-1 


0.03 



This should give a more realistic estimate of decay rates. Table 3.4.4 shows the 
overlap values for a variety of choices of th, hard-core radii, and aB- Also included 
are the results with the MS wavefunction. 
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Empirical Limit on Wavefunction Overlap 

As discussed in section 3.4.3, the H can be lighter than 2 A's without conflicting 
with hypernuclear experiments if it is sufficiently compact, as suggested by some 
models. The constraint imposed by the hypernuclear experiments can be translated 
into an empirical upper limit on the wavefunction overlap between an H and two 
baryons. Using eq. (3.34) for the formation time Tform of an in a double-A 
oxygen- 16 hypernucleus we have 



where /form = $2(mg=2GeV) departure of the phase space factor for hyper- 

nuclear H formation appearing in eq. (3.33), from its value for mu — 2 GeV. 
By crossing symmetry the overlap amplitudes \M.\H^Kk and |A1|AA^_ff only differ 
because the A's in the former are asymptotically plane waves while for the latter 
they are confined to a nucleus; comparing eqns. (3.45) and (3.24) we obtain: 

|A^I^-.AA=^|A^liA^^.- (3.47) 

For oxygen-16, GeV^. Using eqns. (3.46) and (3.47) will give us an 

upper limit on the overlap for the lifetime calculations of the next section. 

Decay rates and lifetimes 

Starting from \M\h-^k\ we can calculate the rates for H decay in various channels, 
as we did for nuclear conversion in the previous section. The rate of if — > decay 
is 



H^nn « K' \' ' ^2{mH) (3.48) 



-2 

2 niH 

X {\M\\^j^\M\h^k\Y , 



47 



where $2 is the phase space factor defined for H ^ nn normahzed as in [55]. Using 
eqs. (3.47) and (3.46), the hfetime for H ^ nn is 

th^nn ~ 9(4) 10^ /^o yr, (3.49) 

for niH — 1-9 (2) GeV, where > 1 \s defined to be (rform/form)/(10^^°s) x 
(5 10^ ^bg)/^- "^^^ ^ therefore cosmologically stable, with a hfetime longer 
than the age of the Universe, if \MWj^_^jj is 10^~^ times smaller than needed to 
satisfy double hypernuclear constraints. As can be seen in Fig. 3.2, this corresponds 
to ^ 1/3 in the IK model discussed above. Note that K1440 and the sensitivity 
to the wavefunction overlap has been eliminated by using rform- 

If niN + rn/s, (2.05 GeV) < niH < 2mA (2.23 GcV), H decay requires only a 
single weak interaction so the rate in eq. (3.48) must be divided by given 
in eqn (3.32). Thus we have 

th-^na ~ 10 iio s. (3.50) 
Finally, if m^j > 2m\ (2.23 GeV), there is no weak interaction suppression and 

TH^AA ^ 4 10-^^ fio s. (3.51) 

Equations (3.49)-(3.51) with //q = 1 give the lower bound on the H lifetime, 
depending on its mass. This result for the H lifetime differs sharply from the 
classic calculation of Donoghue, Golowich, and Holstcin [74], because we rely on 
experiment to put an upper limit on the wavefunction overlap \M.\h-*aa- 
treatment of the color-fiavor-spin and weak interaction parts of the matrix ele- 
ments is approximate, but it should roughly agree with the more detailed calcula- 
tion of ref. [74], so the difference in lifetime predictions indicates that the spatial 
overlap is far larger in their bag model than using the IK and Bethe-Goldstone 
or Miller- Spencer wavefunctions with reasonable parameters consistent with the 
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hypernuclear experiments. The bag model is not a particularly good description of 
sizes of hadrons, and in the treatment of [74] the H size appears to be fixed implic- 
itly to some value which may not be physically realistic. Furthermore, it is hard to 
tell whether their bag model analysis gives a good accounting of the known hard 
core repulsion between nucleons. As our calculation of previous sections shows, 
these are crucial parameters in determining the overlap. The calculation of the 
weak interaction and color-flavor-spin matrix elements in ref. [74] could be com- 
bined with our phenomenological approach to the spatial wavefunction overlap to 
provide a more accurate yet general analysis. We note that due to the small size 
of the H, the p-wave contribution should be negligible. 



3.5 Binding of flavor singlets to nuclei 

After calculating the constraints implied on the H from nuclear transition pro- 
cesses we turn to the calculation of H nuclear binding, ref. [31]. In section 3.3 
we concluded that the relevant experimental constraints on exotic nuclei can place 
strong constraints on the abundance of the H in the case it binds to nuclei. In this 
section we explore binding of flavor singlet to nuclei. We summarize the theory 
of nuclear binding in subsection 3.5.1, to set the framework for and to make clear 
the limitations of our computation. In subsection 3.5.2 we analyze the binding of 
a flavor singlet scalar to nuclei, and calculate the minimum values of coupling con- 
stants needed for binding. Corresponding limits on nucleon-i7 scattering are given 
in subsection 3.5.3. Other flavor-singlets are also considered, in subsection 3.5.4 
and elsewhere. We summarize the results and give conclusions in section 3.6. 
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3.5.1 Nuclear binding-general 

QCD theory has not yet progressed enough to predict the two nucleon inter- 
action ab initio. Models for nuclear binding are, therefore, constructed semi- 
phenomenologically and relay closely on experimental input. 

The long range part of the nucleon-nucleon interaction (for distances r > 1.5 
fm) is well explained by the exchange of pions, and it is given by the one pion 
exchange potential (OPEP). The complete interaction potential Vy is given by 

~'~ ^ij^ where vfj contains all the other (heavy meson, multiple meson and quark 
exchange) parts. In the one boson exchange (OBE) models the potential vfj arises 
from the following contributions: 

• In the intermediate region (at distances around r ~ 1 fm) the repulsive 
vector meson (p, uj) exchanges are important. A scalar meson denoted a was 
introduced to provide an attractive potential needed to cancel the repulsion 
coming from the dominant vector cu meson exchange in this region. Moreover, 
a spin-orbit part to the potential from both a and cu exchange is necessary 
to account for the splitting of the phase shifts in NN scattering. 

• At shorter scales (r < 1 fm) , the potential is dominated by the repulsive 

vector meson (p, cu) exchanges. 

• For r ^ 0.5 fm a phenomenological hard core repulsion is introduced. 

However, many of these OBE models required unrealistic values for the meson- 
nucleon coupling constants and meson masses. With this limitation the OBE the- 
ory predicts the properties of the deuteron and of two-nucleon scattering, although, 

it cannot reproduce the data with high accuracy. 

A much better fit to the data is obtained by using phenomenological potentials. 
In the early 1990's the Nijmegen group [75] extracted data on elastic NN scattering 
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and showed that all NN scattering phase shifts and mixing parameters could 
be determined quite accurately. A^A^ interaction models which fit the Nijmegen 
database with a x^/^data ~ 1 are called 'modern'. They include Nijmegen models 
[76], the Argonne Vig [77] and CD-Bonn [78] potentials. These potentials have 
several tens of adjustable parameters, and give precision fits to a wide range of 
nucleon scattering data. 

The construction of 'modern' potentials can be illustrated with the Nijmegen 
potential. That is an OBE model based on Regge pole theory, with additional 
contributions to the potential from the exchange of a Pomeron and f, f ' and A2 
trajectories. These new contributions give an appreciable repulsion in the central 
region, playing a role analogous to the soft or hard core repulsion needed in semi- 
phenomenological and OBE models. 

Much less data exists on hyperon-nucleon interactions than on NN interactions, 
and therefore those models are less constrained. For example the extension of the 
Nijmegen potential to the hypcr-nuclear (YN) sector [79] leads to under-binding 
for heavier systems. The extension to the AA and 'EN channels cannot be done 
without the introduction of extra free parameters, and there are no scattering data 
at present for their determination. 

The brief review above shows that the description of baryon binding is a difficult 
and subtle problem in QCD. Detailed experimental data were needed in order to 
construct models which can describe observed binding. In the absence of such 
input data for the H analysis, we must use a simple model based on scalar meson 
exchange described by the Yukawa potential, neglecting spin effects in the nucleon 
vertex in the first approximation. We know from the inadequacy of this approach 
in the NN system that it can only be used as a crude guide. However since the 
strength of couplings which would be needed for the H to bind to light nuclei are 
very large, compared to their expected values, we conclude that binding is unlikely. 
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Thus limits on exotic nuclei cannot be used to exclude the existence of an H or 
other compact flavor singlet scalar or spin- 1/2 hadron. 

3.5.2 Binding of a flavor singlet to nuclei 

The H cannot bind through one pion exchange because of parity and also fla- 
vor conservation. The absorption of a pion by the H would lead to an isospin 
I — 1 state with parity (— 1)"^+^, which could be AE'^ or heavier Sp compos- 
ite states. These states have mass ^0.1 GeV higher than the mass of the H 
(for rriA + niN ^'^Tn\), which introduces a strong suppression in 2"^ order 

perturbation theory. Moreover, the baryons in the intermediate state must have 
relative angular momentum L = 1, in order to have odd parity as required; this 
introduces an additional suppression. Finally, production of AE° or EN states is 
further suppressed due to the small size of the H, as explained in §3.4. Due to all 
these effects, we conclude that the contribution of one or two pion exchange to H 
binding is negligible. 

The first order process can proceed only through the exchange of a fiavor singlet 
scalar meson and a glueball. The hghtest scalar meson is f (400-1 100) (also called 
a). The mass of the glueball is considered to be around ~ 1.5 GeV. In Born 
approximation, the Yukawa interaction leads to an attractive Yukawa potential 
between nucleons 

V(r) = -^-e-''^ (3.52) 
An r 

where ii is the mass of the exchanged singlet boson s (a or glueball) and gg' is the 
product of the s-H and s-nucleon coupling constants, respectively. The potential of 
the interaction of if at a position r with a nucleus, assuming a uniform distribution 
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of nucleon p = 4 inside a nuclear radius R, is then 



gg'A 



47r V 



-^\r—r'\ 



v^-^-:- / ^ — ^rfV, (3.53) 



where A is the number of nucleons, V is the volume of the nucleus and r is the 
position vector of the H. After integration over the angles the potential is 



where we used R — 1.3574^/^ fm; 



2/1 



1 - (1 + i^R) e-'^^ 2i2MH 



r < R 



2// [iiRcosh[iiR] - sinh[//i?]] ^ r > R. 



Throughout, we use h—c—1 when convenient. 

Fig. 3.4 shows the potential the nucleus presents to the H for A = 50, taking the 
mass of the exchanged boson to be = 0.6 and 1.5 GeV. The depth of the potential 
is practically independent of the number of nucleons and becomes shallower with 
increasing scalar boson mass /i. 

Note that Born approximation is applicable at low energies and for small cou- 
pling constants; it may not be valid for H binding. Born approximation is valid 
when 

m qq' 

-f- « 1, 3.55 

where m is the reduced mass and the mass of the exchanged particle. As we shall 
see, this condition is actually not satisfied for values of gg' which assure binding for 
the i^-mass range of interest. This underlines the fact that no good first principle 
approach to nuclear binding is available at present. 

We can now calculate the value of c* = j for which the potential is equal 
to the minimum value needed for binding; in square well approximation this is 
given by 

Vmin = ^-5^- (3-56) 
8R^m 
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Figure 3.4: Potential in GeV, for ^=1, A=50 and /i = 0.6 (dashed) or = 1.5 
GeV (solid) as a function of distance r. 

Fig. 3.5 shows the dependence of c* on the mass of the exchanged particle, /x. The 
maximum value of for which the H does not bind decreases with increasing H 
mass, and it gets higher with increasing mass of the exchanged particle, /i. 

The H does not bind to hght nuclei with ^4 < 4, as long as the product of 
couplings 

< [0.27, 0.73, 1.65], for ^ = [0.6, 1, 1.5] GeV, (3.57) 

where c = gNNa gHHc/{^'^) or Qnng gHHo/i^'^)- The H will not bind to heavier 
nuclei if 

< [0.019, 0.054, 0.12], for = [0.6, 1, 1.5] GeV. (3.58) 

In the next sections we will compare these values to expectations and limits. 

It should also be noted that binding requires the product of coupling constants, 
gg' to be positive and this may not be the case. Even in the case of hyperons, 

experimental information was necessary to decide whether the S has a relative 
positive coupling [80]. 
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Figure 3.5: Critical value c* of the coupling constant product versus nuclear size 
needed for the H to just bind, for /i[GeV]= 0.7 (dotted), 1.3 (dashed) and 1.5 
(solid). 

3.5.3 Limits on cm from Nucleon H elastic scattering 

The nucleon- if elastic scattering cross section is expected to be very small, due 
to the compact size of the H and the suppression of color fluctuations on scales 
^ 1 GeV~^ in the nucleon. Ref. [40] estimates ~ lO^""* mb. This can be trans- 
lated to an estimated upper limit on the product c m which determines the po- 
tential presented to the if by a nucleus, as follows. In the one boson exchange 
model, the elastic H-N cross section due to the a- or glueball-mediated Yukawa 
interaction is given by 



da 



{2mcf 




1 



(3.59) 



In the low energy limit 




(3.60) 



Writing (Thn — C-3IO ^ mb and = A*Gevl GeV, this gives 



cm — 0.007^(7_3 //q, 



GeV. 



(3.61) 
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Comparing to the values of c* needed to bind, we see that for mn < 2mp this is 
too small for the H to bind, even to heavy nuclei^. 

If dark matter consists of rehc if's, we can demonstrate that if's do not bind 
to nuclei without relying on the theoretical estimate above for anN- It was shown 
in [81] that the XQC experiment excludes a dark matter-nucleon cross section 
axN larger than about 0.03 mb for mx ~ 1-5 GeV. Thus if dark matter consists 
of a stable H it would require axN < 0.03 mb, implying c < [0.01,0.03,0.06] for 
fjL — [0.6, 1.0, 1.5] GeV and the H would not bind even to heavy nuclei. 

A generic new scalar flavor singlet hadron X which might appear in an extension 
of the Standard Model, might not have a small size and correspondingly small value 
of (JxN, and it might not be dark matter and subject to the XQC limit. In that 
case, it is more useful to turn the argument here around to give the maximum 
above which the X would bind to nuclei in the OBE approximation. Prom eqn 
(3.54), (3.5.2) and /(O) = 2// we have 

_ 7r^(1.35 fm),.^ 

Then eq. (3.60) leads to 

a*xj^ ^ 155 ^-^/^ mb. (3.63) 

That is, for instance, if it is required that the X not bind to He then it must have 
a cross section on nucleons lower than 25 mb. 



3.5.4 Flavor singlet fermion 

The analysis of the previous sections can be extended to the case of a flavor singlet 
fermion such as the glueballino - the supersymmetric partner of the glueball which 
^We have summarized the net effect of possibly more than one exchange boson {e.g., a and 
glueball) by a single effective boson represented by a cf^ and /Xeff • 
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appears in theories with a hght gluino [82]. In this case the possible exchanged 
bosons includes, in addition to the cr and the glueball, the flavor-singlet component 
of the pseudoscalar meson 77' (m^ = 958 MeV). However the size of the Rq should 
be comparable to the size of the glueball, which was the basis for estimating the 
size of the H. That is, we expect rji^ ^ ^ and thus aji^^ ~ o'hn [40]. Then 
arguments of the previous section go through directly and show the Rq is unlikely 
to bind to light nuclei^. 

3.6 The Existence of the H - Summary 

In §3.3 we considered the constraints placed on the H dibaryon by the stability of 
nuclei and hypcrmiclei with respect to conversion to an H, and we calculated the 
lifetime of the H in the case when it is heavier than two nucleons. In the model 
calculations we used the Isgur-Karl wavefunctions for quark in baryons and the H, 
and the Miller-Spencer and Bruecker-Bethe-Goldstone wavefunctions for nucleons 
in a nucleus, to obtain a rough estimate of the /J-baryon-baryon wavefunction over- 
lap. By varying the IK oscillator strength parameter and the hard-core radius in 
the BBG wavefunction over extreme ranges, we find that the wavefunction overlap 
is very sensitive to the size and shape of the hadronic and nuclear wavefunctions. 
With the BBG (MS) wavefunction, the hypernuclear formation time of an H is 
comparable to or larger than the decay time for the A's and thus the H is not 

excluded, if TH^lf^ (1/3) r^. ^ We conclude that the observation of A decays 

■'Niissinov [83] considered that the Rq would bind to nuclei, by assuming that the depth of 

the potential presented by the nucleus to the Rq is at least 2-4 MeV for 16 < A < 56. However 

the discussion of the previous sections, with (Jr^n = 10~^cr_3 mb, gives a potential depth of 0.07 

MeV V^/(™flo/GcV). 

^Thc overlap between an H and two nucleons should be strongly suppressed also in the Skyrme 

model, in view of the totally different nature of the nucleon and H solitons [84,85]. However, a 
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in double- A hypernuclei cannot be used to exclude mn < 2mA, given our present 
lack of understanding of the hadronic and nuclear wavef unctions. 

In the second part of our work we abstracted empirical relations which give us 
relatively model-independent predictions for the H lifetime. By crossing symmetry, 
the overlap of the wavefunctions of an H and two baryons can be constrained using 
experimental limits on the formation time of an if in a hypernucleus. Using the 
empirically constrained wavefunction overlap and phenomenologically determined 
weak interaction matrix elements, we can estimate the lifetime of the H with 
relatively little model uncertainty. We find: 

(i) If niH > 2mA, the hypernuclear constraint is not applicable but the H would 
still be expected to be long-lived, in spite of decaying through the strong 
interactions. E.g., with the BBG wavefunction and ri/ ^ 1/2 tat, tj^ > 4 10~^^ 
sec. 

(ii) If mAT -l- mA ^ mn ^ 2mA, the H lifetime is ^ 10 sec. 

(iii) If 2mN^rnH^rnN + rriA, the H lifetime is > 10^ yr. For (1/3) as 
suggested by some models, the H lifetime is comparable to or greater than 
the age of the Universe. 

Our results have implications for several experimental programs: 

(i) The observation of A decays from double A hypernuclei excludes that Tform, 
the formation time of the if in a double A hypernucleus, is much less than 
ta. However if rform is of order ta, some double A hypernuclei would produce 
an H. One might hope these if's could be observed by reconstructing them 
through their decay products, e.g., H E~p. Unfortunately, our calculation 
method for computing the overlap has not been developed so we are unable to explore this here. 
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shows that th ^ 10 sec for the relevant range of mn, so any if's produced 
would diffuse out of the apparatus before decaying. 

(ii) Some calculations have found rriH < 2{mp + rrie), in which case the H is ab- 
solutely stable and nucleons in nuclei may convert to an i?. We showed that 
SuperK can place important constraints on the conjecture of an absolutely 
stable H, or conceivably discover evidence of its existence, through observa- 
tion of the pion(s), positron, or photon produced when two nucleons in an 
oxygen nucleus convert to an H. We estimate that SuperK could achieve a 
lifetime limit r ^ few 10^^ yr. This is the lifetime range estimated with the 
BBG wavefunction for mij> 1740 MeV and rn ~ 1/5 r^v- An H smaller 
than this seems unlikely, so mn ^ 1740 MeV is probably already ruled out. 

In §3.5.1 we first reviewed the theory of nuclear binding and emphasized that 
even for ordinary nucleons and hyperons there is not a satisfactory first-principles 
treatment of nuclear binding. We showed that exchange of any pseudoscalar meson, 
or of two pseudoscalar octet mesons, or any member of the vector meson octet, 
makes a negligible contribution to the binding of an H or other flavor singlet 
scalar hadron to a nucleon. The dominant attractive force comes from exchange 
of a glueball or a cr (also known as the f (400-1 100) meson), which we treated with 
a simple one boson exchange model. The couplings of a and glueball to the H 
are strongly constrained by limits on anN, to such low values that the H cannot 
be expected to bind, even to heavy nuclei. Thus we conclude that the strong 
experimental limits on the existence of exotic isotopes of He and other nuclei do 
not exclude a stable H. 

More generally, our result can be applied to any new flavor singlet scalar particle 
X, another example being the supersymmetric hybrid baryon (udsg) discussed 
in [82]. If axN < 25 mb GeV/mx, the X particle will not bind to light nuclei 
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and is "safe". Conversely, if ax at >> 25 mb GeV/mx, the X particle could 
bind to light nuclei and is therefore excluded unless if there is some mechanism 
suppressing its abundance on Earth, or it could be shown to have an intrinsically 
repulsive interaction with nucleons. This means the self-interacting dark matter 
(SIDM) particle postulated by Spcrgel and Steinhardt [86] to amehorate some 
difficulties with Cold Dark Matter, probably cannot be a hadron. SIDM requires 
(^xx/Mx 0.1 — 1 b/GeV; if X were a hadron with such a large cross section, then 
on geometric grounds one would expect uxn ~ l/4crxx which would imply the X 
binds to nuclei and would therefore be excluded by experimental limits discussed 
above. 



As we have seen in the previous sections, the existence of the H di-baryon is not 
ruled out by experiment if the H is tightly bound and compact. If mn ^ttin + '>tia 
and r//^l/3 the H can be cosmologically stable. In the rest of this section 
we explore the possibility that the DM consists of the H or H,H and therefore be 
predicted within the Standard Model. 

1. The H Dark Matter. The number density of nonrelativistic species i in 
thermal equilibrium is given by 



If baryon number is conserved we have ij,h = 2/X7v, and therefore the nucleon 
and the H energy densities are related as 



3.7 The H or H, H as Dark Matter 




(3.64) 



mnriH 



(27r) 
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2m„ - ruH 



T 



(3.65) 
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The left-hand side is actually the ratio Qh/^n and for H DM it is fixed, see 
eq. (2.1). Then, by solving eq. (3.65) we can find Tf.o., the temperature at 
which the H has to go out of equilibrium in order to have the proper DM 
energy density today. The equation has a solution only for uih ^'^ttin, a 
mass range which is unfavored by the discussion in §3.4.3. The freeze-out 
temperatures are 15 (7) MeV, for rriH = 1-5 (1.7) GcV. These temperatures 
correspond to an age of the Universe of 0.003 (0.02) sec. By that time all 
strangeness carrying particles already decayed and therefore the H cannot 
stay in equilibrium (for example, through reactions as K^H <-> pA) at such 
low temperatures. We see that even if the H was cosmologically stable it 
could not be thermally produced in the early universe in sufficient amount 
to be dark matter. 

2. The HH Dark Matter. In this case the H would decouple as in the B- 
sequestration scenario discussed in §2. The reactions which would keep it in 
the equilibrium, up to the temperatures listed in Table 2.1, are 

HN ^ KK 
HN ^ KK. 

In this scenario As and Ks stay in equilibrium sufficiently long, because the 
above reactions proceed in the left direction with the rates which are, for 
temperatures of interest, much higher than the decay rates of A and K. The 
H and H stay in equilibrium through these interactions and may reach DM 
abundance. 

In the next Chapter we explore direct and indirect experimental constraints on the 
HH DM scenario. 
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Chapter 4 

Dark Matter constraints in the 
range of the H parameters 

This Chapter is organized as follows: we calculate direct DM detection experiments 
in the mass and cross section range expected for the i7 in § 4.1 and indirect 
constraints that could be place on H DM in § 4.2; we show that the H DM could 
be ruled out from the heat production in Uranus. 

4.1 Direct DM Detection - Light DM Constraints 

In this Section we focus on direct detection experiments in order to see whether 
the H is allowed as a DM candidate given that his expected cross section with 
nucleons is of the order //b (see discussion in § 3.2). More generally we explore the 
limits on light DM, with primary interest in //b cross section range. 

The region of mass m > 10 GeV is well explored today, up to TeV range, from 
strong (cT ~ 10 mb) to weak (a ~ 10~^pb) cross sections. Here we explore the 
possibility that the DM mass is in the range 0.4<m<10 GeV. Masses below 



62 



~ 0.4 GeV are below the threshold of direct detection DM experiments and are 
therefore unconstrained, with the exception of axion searches. 

The mass range m < 10 GeV has not yet been explored carefully. Several dark 
matter underground experiments have sufficiently low mass threshold today: the 
CRESST [87], DAMA [88], IGEX [89], COSME [90] and ELEGANT [91] experi- 
ments. Except for DAMA, these experiments have published upper limits on the 
cross section assuming it is weak, but have not addressed the case of stronger cross 
sections,^ where the approach for extracting the cross section limits is substan- 
tially different, as we explain below. Also, recent data from an X-ray experiment 
XQC, [93] proved to be valuable in constraining low mass DM, but limits based on 
the final data have not yet been derived. Since XQC is a space- based experiment 
it is especially suitable for exploring the higher cross section range. In [81] it was 
shown that in the low mass range the XQC experiment rules out Strongly Inter- 
acting DM (SIMPs, [86]). Dark matter with low masses and 'intermediate' cross 
sections, several orders of magnitude smaller than normal hadronic cross sections, 
remains to be fully analyzed and that is the focus of this work. We will abbreviate 
DM with intermediate cross section on nucleons as DMIC. 

Early limits from DMIC direct detection experiments can be found in the pa- 
per [94] by Rich, Rocchia and Spiro in which they reported results from a 1987 
balloon experiment. Starkman et al. [95] reviewed DM constraints down to a 
mass of 1 GcV as of 1990. Wandelt et al. [81] added constraints based on pre- 
liminary data from the more recent XQC sounding rocket experiment. The above 
constraints are discussed and updated further in the text. In previous works on 
this topic the effect of the halo particle's velocity distribution on the cross section 
limit was not explored. Since the only detectable light particles arc those in the 

"'^DAMA did publish strong cross section limits in [92], but they were based on a dedicated 
experiment which had a higher mass threshold to > 8 GeV. 
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exponential tail of the velocity distribution, the limits on light DM are sensitive 
to the parameters in the velocity distribution, in particular to the value of the 
escape velocity cutoff. We investigate this sensitivity in the standard Isothermal 
Sphere model, where the DM velocity distribution function is given by a truncated 
Maxwell-Boltzmann distribution. We also consider the spin-independent and spin- 
dependent interaction cases separately. Except in Section 4.1.5, we assume a single 
type of DM particle. 

4.1.1 Direct Dark Matter Detection 

The basic principle of DM detection in underground experiments is to measure the 
nuclear recoil in elastic collisions, see for example [96]. The interaction of a DM 
particle of mass m < 10 GeV, produces a recoil of a nucleus of 20 keV or less. The 
recoil energy (which grows with DM mass as in (2) below) can be measured using 
various techniques. The detectors used so far are: 

• ionization detectors of Ge (IGEX, COSME) and of Si 

• scintillation crystals of NaT (DAMA, ELEGANT) 

• scintillators of CaF2 

• liquid or liguid gas Xenon detectors (DAMA, UKDMC) 

• thermal detectors (bolometers) with saphire (CRESST, ROSEBUD) , telurite 
or Ge (ROSEBUD) absorbers 

• bolometers, which also measure the ionization like that of Si (CDMS) and 
Ge (CDMS, EDELWEISS). 

As we have seen in the Introduction, the most accepted DM candidate in GeV 
range is the neutralino. The expected signal for neutralino-nucleon scattering is in 
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the range 10 to 10^^ counts/kgday, [97]. The smallness of signal dictates the ex- 
perimental strategies, detectors have to be highly radio pure, and have substantial 
shielding (they are built deep underground). 

For a given velocity distribution f{v}, the differential rate per unit recoil energy 
Eji in (kg day keV)^^ in the detector can be expressed as 



where nx is the number density of DM particles, Nt is the number of target nuclei 
per kg of target, axA is the energy dependent scattering cross section of DM on a 
nucleus with mass number A, g{v) is the probability that a particle with velocity 
V deposits an energy above the threshold Eth in the detector, and Vmin is the 
minimum speed the DM particle can have and produce an energy deposit above 
the threshold. The recoil energy of the nucleus is given by 



where 6cm is the scattering angle in the DM-nucleus center of mass frame. We 
will assume isotropic scattering as is expected at low energies. So, for instance, for 
A — 16, m — 1 GeV and an energy threshold of 600 eV, the minimal DM velocity 
to produce a detectable recoil is Vmin = 680 km/s, in the extreme tail of the DM 
velocity distribution. 

In order to compare cross section limits from different targets we will normalize 
them to the proton-DM cross section, axp- For the simplest case of interactions 
which are independent of spin and the same for protons and neutrons, the low 
energy scattering amplitude from a nucleus with mass number A is a coherent sum 
of A single nucleon scattering amplitudes. The matrix element squared therefore 
scales with size of nucleus as ~ In addition the kinematic factor in the cross 




(4.1) 




(4.2) 
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section depends on the mass of the participants in such a way [96, 98] that 



SI r / .Nn2 



A' (4.3) 



where fi{A) is the reduced mass of the DM-nucleus system, and fi{p) is the reduced 
mass for the proton-DM system. At higher momentum transfer = 2mjv-£'R the 
scattering amphtudes no longer add in phase, and the total cross section (Txa(?) 
becomes smaller proportionally to the form factor F'^{q^), axAio) — Co-^^(?^)- 

We take this change in the cross section into account when we deal with higher 
mass (m ^ 10 GcV) dark matter; for smaller masses the effect is negligible. We 
adopt the form factor F{q^) = exp [— l/10(gi?)^] with R = 1.2^4^/^ fm, used also 
in [99,100]. The simple exponential function is suffitiently accurate for our pur- 
poses and easy to implement using the Monte Carlo method to sample momentum 
transfer q, from its distribution given by the form factor. The procedure is de- 
scribed in more detail in Appendix A. 

For spin dependent interactions the scattering amplitude changes sign with the 
spin orientation. Paired nucleons therefore contribute zero to the scattering ampli- 
tude and only nuclei with unpaired nucleon spins are sensitive to spin dependent 
interactions. Due to the effect of coherence, the spin independent interaction is usu- 
ally dominant, depending on the mass of the exchanged particle [101]. Therefore, 
the spin dependent cross section limit is of interest mainly if the spin independent 
interaction is missing, as is the case, for example, with massive majorana neutri- 
nos. Another example of DM with such properties is photino dark matter, see [98], 
in the case when there is no mixing of left- and right- handed scalar quarks. The 
amplitude for DM-nucleus spin dependent interaction in the case of spin 1/2 DM, 
in the nonrelativistic limit, is proportional to [98, 102] 

M^{N\J\N)-sx (4.4) 
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where J is the total angular momentum of the nucleus, \N > are nuclear states 
and sx is the spin of the DM particle. In the case of scalar DM the amplitude is 

M-^ {N\J\N)-{qx^) (4.5) 

where q and ^ are the initial and final momenta of the scattering DM particle. 
Thus the cross section for this interaction is proportional to the fourth power of 
the ratio q/M, of DM momentum to the mass of the target which enters through 
the normalization of the wavefunction. Therefore the spin dependent part of the 
interaction for scalar DM is negligible when compared to the spin independent 
part. 

We adopt the standard spin-dependent cross section parametrization [96] 

TxA ~ ^iAf [X\J{J + 1)UC],^ (4.6) 

where A is a parameter proportional to the spin, orbital and total angular momenta 
of the unpaired nucleon. The factor C is related to the quark spin content of the 
nucleon, C = Xl^l^g' ^ ~ u,d,s, where T^''''* is the charge of the quark type q 
and Ag is the fraction of nucleon spin contributed by quark species q. The nuclear 
cross section normalized to the nucleon cross section is 

' [AV(J + 1)U (CxaV .... 
[AV(J + 1)], \CxJ ' ^ ■ ^ 

The values of proton and neutron C factors, Cxp, Cxn vary substantially depend- 
ing on the model. For targets of the same type - odd-n (Si, Ge) or odd-p (Al, Na, 
I) nuclei - this model dependence conveniently cancels. The comparison of cross 
sections with targets of different types involves the Cxp/Cxn ratio. This ratio was 
thought to have the value ~ 2 for any ncutralino, based on the older European 
Muon Collaboration (EMC) measurements, but the new EMC results imply a ratio 
which is close to one for pure higgsino, and is ^ 10 otherwise. (The biggest value 
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for the ratio is Cp/Cn ~ 500, for bino.) We normalize our spin dependent results 
to the proton cross section uxp using Cxp/Cxn = 1 for definiteness below. 

In this paper we assume that the DM halo velocity distribution is given by a 
truncated Maxwell- Boltzmann distribution in the galactic reference frame, as in 
the Isothermal Sphere Halo model [103]. We direct the z axis of the Earth's frame 
in the direction of the Local Standard of Rest (LSR) motion. ^ The DM velocity 
distribution, in the Earth's frame, is given by 



f{vz,v±) = Nexp 



(4.8) 



Here Vc is the local circular velocity and it is equal to \/2 times the radial velocity 
dispersion in the isothermal sphere model; ve is the velocity of the Earth in the 
Galactic reference frame. Throughout, superscript "t" indicates a tangential com- 
ponent. This neglects the Earth's motion in the radial direction which is small. 
The velocities Vz and v\_ are truncated according to \Jv1 + v\_ ^Vesc, where v^sc is 
the escape velocity discussed below. 

The model above is the simplest and the most commonly used model which 
describes a self-gravitating gas of coUisionless particles in thermal equilibrium. On 
the other hand numerical simulations produce galaxy halos which are triaxial and 
anisotropic and may also be locally clumped depending on the particular merger 
history (see [104] for a review). This indicates that the standard spherical isotropic 
model may not be a good approximation to the local halo distribution. Here we 
aim to extract the allowed DM window using the simplest halo model, but with 
attention to the sensitivity of the limit to poorly determined parameters of the 
model. The effects of the more detailed halo shape may be explored in a further 
work. 



^The Local Standard of Rest used here is the dynamical LSR, which is a system moving in a 
circular orbit around the center of Milky Way Galaxy at the Sun's distance. 
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We ignore here the difference between the DM velocity distribution on the 
Earth, deep in the potential well of the solar system, and the DM velocity distri- 
bution in free space. This is a common assumption justified by Gould in [105] as 
a consequence of Liouville's theorem. Recently Edsjo et al. [29] showed that the 
reahstic DM velocity distribution differs from the free space distribution, but only 
for velocities km/s. Therefore, the free space distribution is a good approx- 

imation for our analysis, since for light dark matter the dominant contribution to 
the signal comes from high velocity part of the distribution. 

The velocity of the Earth in the Galactic reference frame is given by 

Ve = VlSR + Vs + VE,orb: (4-9) 

where vlsr is the velocity of the local standard of rest LSR: it moves with local 
circular speed in tangential direction v^^^ = Vc, toward I — 90°, b — 0°, where 
/ and b are galactic longitude and latitude. The velocity of the Sun with respect 
to the LSR is vs = 16.6 km/s and its direction is / = 53", b = 25° in galactic 
coordinates. VE,orb — 30 km/s is the maximal velocity of the Earth on its orbit 
around the Sun. 

The magnitude of vl^gj^ has a considerable uncertainty. We adopt the conserva- 
tive range Vc = (220±50) km/s which relies on purely dynamical observations [106]. 
Measurements based on the proper motion of nearby stars give a similar central 
value with smaller error bars, for example Vc{Ro) — (218±15) km/s, from Cepheids 
and Vc{Ro) — (241 ± 17) km/s, from SgrA* (see [107] and references therein). The 
choice Vc — (220 ± 50) km/s is consistent with the DAMA group analysis in [108] 
where they extracted the dependence of their cross section hmits on the uncertainty 
in the Maxwellian velocity distribution parameters. 



69 



Projecting the Earth's velocity on the tangential direction (/ = 90°, b = 0°) we 

get 

vl;^Vc + vl + v\^^^ co%\u{t - to)] (4.10) 

where v\ = 12 km/s; = 30 cos 7 km/s where cos 7 = 1/2 is the cosine of 
the angle of the inclination of the plane of the ecliptic, u = 27r/365 day^^ and 
^0 is June 2nd, the day in the year when the velocity of the Earth is the highest 
along the LSR direction of motion. In the course of the year cos[uj{t — to)] changes 
between ± 1, and the orbital velocity of the Earth ranges ±15 km/s. Taking all of 
the uncertainties and annual variations into account, the tangential velocity of the 
Earth with respect to the Galactic center falls in the range = (167 to 307) km/s. 

The other parameter in the velocity distribution with high uncertainty is the 
escape velocity, Vesc — (450 to 650) km/s [109]. We will do our analysis with the 
standard choice of velocity distribution function parameters, 

v% = 230 km/s, Vc = 220 km/s, Vesc = 650 km/s, (4.11) 

and with the values oi ve and Vgsc from their allowed range, which give the lowest 
count in the detector and are therefore most conservative: 

= 170 km/s, = 170 km/s, Ve.c = 450 km/s. (4.12) 

For experiments performed in a short time interval we take the vahic of ^,^f^ cos[c<j(t— 
to)] which corresponds to the particular date of the experiment, and the lowest 
value of v% allowed by the uncertainties in the value of Vc- 

Another effect is important in the detection of moderately interacting particles. 
Since particles loose energy in the crust rapidly (mean free path is of the order of 
100 m) only those particles which come to the detector from 27r solid angle above 
it can reach the detector with sufficient energy. Since the velocity distribution of 
the particles arriving to the detector from above depends on the detector's position 
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on Earth with respect to the direction of LSR motion, the detected rate for a hght 
DMIC particle will vary with the daily change of position of the detector. This 
can be a powerful signal. 

4.1.2 XQC Experiment 

For hght, moderately interacting dark matter the XQC experiment places the 
most stringent constraints in the higher cross section range. The XQC experiment 
was designed for high spectral resolution observation of diffuse X-ray background 
in the 60 — 1000 eV range. The Si detector consisted of an array of 36 1 mm^ 
microcalorimeters. Each microcalorimeter had a 7000 Angstrom layer of HgTe X- 
ray absorber. Both the HgTe and the Si layers were sensitive to the detection. The 
experiment was performed in a 100 s flight in March, and therefore the Earth's 
velocity falls in the 200 to 300 km/s range. The experiment was sensitive 
to energy deposit in the energy range 25 — 1000 eV. For energy deposits below 
25 eV the efficiency of the detector drops off rapidly. For energy deposits above 
about 70 eV the background of X-rays increases, so XQC adopted the range 25-60 
eV for extraction of DM hmits, and we will do the same. This translates into a 
conservative mass threshold for the XQC experiment of 0.4 GeV, obtained with 
Vesc = 450 km/s and = 200 km/s, which is the lowest mass explored by direct 
DM detection apart from axion searches. 

The relationship between the number of signal events in the detector Ns and 
the scattering cross section axA of DM particles on nuclei is the following 

Ns^nx fT (Arsi(i;si)(7si + A^Hg[(^'Hg)(7Hg + {vTeWTe)], (4.13) 

where N^i and A^Hg are the numbers of Si and Hg (Te) nuclei in the detector, 
nx is the local number density of DM particles, (vsi), {vug) and {vTe) are the 
effective mean velocities of the DM particles on the Si and HgTe targets, / is 
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the efficiency of the detector, and T = 100 s is the data-taking time. In this 
energy range, / ~ 0.5. The standard value for the local DM energy density is 
px — 0.3 GeV cm~^. However, numerical simulations combined with observa- 
tional constraints [110] indicate that the local DM energy density px may have 
a lower value, 0.18<px/(GeV cm~^) ^0.3. In our calculations we use both the 
standard value px = 0.3 GeV/cm^, and the lower value suggested by the numerical 
simulations, px — 0.2 GeV/cm^. The cross sections crsi, CTRg, cxe are calculated 
using equations (4.3) and (4.7). In this section and the next we assume that DM 
has dominantly spin-independent cross section with ordinary matter. In § 4.1.4 
we consider the case of DM which has exclusively spin-dependent cross section or 
when both types of interaction are present with comparable strength. 

XQG observed two events with energy deposit in the 25-60 eV range, and ex- 
pected a background of 1.3 events. The equivalent 90% cl upper limit on the 
number of signal events is therefore Ns — 4.61. This is obtained by interpolating 
between 4.91 for expected background = 1.0 event and 4.41 for expected back- 
ground = 1.5 events, for 2 observed events using table IV in ref. [111]. 

We extract the cross section limits using our simulation. Because of the screen- 
ing by the Earth we consider only particles coming from the 27r solid angle above 
the detector, for which {n- v) < and for them we simulate the interaction in the 
detector, for particles distributed according to (4.8). We take the direction of the 
LSR motion, h as the z axis. 

We choose the nucleus i which the generated DM particle scatters from, using 
the relative probability for scattering from nucleus of type i, derived in Appendix 
B: 

P. = ^ = ^?^^, (4.14) 
where Aj is the mean free path in a medium consisting of material with a mass 
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number Af. Aj = {niaxAi)~^ ■ Here Ui is the number density of target nuclei i in 
the crust, cxAi is the scattering cross section of X on nucleus Aj and the effective 
mean free path, Xeff, is given as 



In each scattering the DM particle loses energy according to (4.2), and we assume 
isotropic scattering in the cm. frame. 

We determine the effective DM velocity < va > as 



where the sum is over the velocities of those DM particles which deposit energy in 
the range 25-60 cV, in a collision with a nucleus of type A, and Ntot is the total 
number of generated DM particles. The result depends on the angle between the 
experimental look direction, and the motion of the Earth. The zenith direction 
above the place where the rocket was launched, uxqc, is toward h — +82°, I — 75°. 
Thus the detector position angle compared to the direction of motion of the Earth 
through the Galaxy is 82°. Only about 10% of the collisions have an energy deposit 
in the correct range. Putting this together gives the 90% confidence level XQC 
upper limit on the spin independent cross section for DMIC shown in Figures 
4.3 and 4.4. The solid hne limit is obtained using the most conservative set of 
parameters (p = 0.2 GeV/cm^, = 200 km/s, Vesc = 450 km/s) and the dotted 
line is the limit obtained by using the standard parameter values in cq. (4.11). The 
upper boundary of the upper domain, a ~ 10^ -j- 10^ mb is taken from [81]. 

When the dark matter mass is higher than 10 GeV, the form factor suppression 
of cross section is taken into account. We give the details of that calculation in 
the Appendix B. 

In the next section we explain how the upper boundaries of the excluded region 
from the underground experiments shown in these figures are obtained. Also shown 




(4.15) 




(4.16) 
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are the original limits from the balloon experiment of Rich et al. [94] obtained using 
the "standard" choices for DM at the time of publishing (dashed line) as well as 
the limits obtained using the conservative values of parameters {v% — 170 km/s, 
since the experiment was performed in October, and Vesc — 450 km/s). Fig. 4.4 
zooms in on the allowed window in the m < 2.4 GeV range. 

4.1.3 Underground Detection 

In this section we describe the derivation of the lower boundary of the DMIC 
window from the underground experiments. This is the value of axp above which 
the DM particles would essentially be stopped in the Earth's crust before reaching 
the detector. (More precisely, they would loose energy in the interactions in the 
crust and fall below the threshold of a given experiment.) To extract the limit 
on axp we generate particles with the halo velocity distribution and then follow 
their propagation through the Earth's crust to the detector. We simulate the 
DM particle interactions in the detector and calculate the rate of the detector's 
response. We compare it to the measured rate and extract cross section hmits. The 
basic input parameters of our calculation are the composition of the target, the 
depth of the detector and the energy threshold of the experiment. We also show 
the value of the dark matter mass threshold mxH, calculated for the standard and 
conservative parameter values given in eq. (4.11) and eq. (4.12). The parameters 
are summarized in Table 4.1.3 for the relevant experiments. 

In the code, after generating particles we propagate those which satisfy (fi-v) < 
through the crust. Given the effective mean free path in the crust eq. (4.15), the 
distance traveled between two collisions in a given medium is simulated as 

X ^ -Xeff In R (4.17) 

where R is a uniformly distributed random number, in the range (0,1). After 
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Table 4.1: The parameters of the experiments used for the extraction of cross 
section hmits; uith is the minimum mass of DM particle which can produce a 
detectable recoil, for the standard and conservative parameter choice. The energy 
threshold values refer to the nuclear response threshold. This corresponds to 
the electron response threshold divided by the quenching factor of the target. 



Experiment 


Target 


Depth 






CRESST, [87] 


AI2O3 


1400 m 


600 eV 


0.8 (1.1) GeV 


DAMA. [88] 


Nal 


1 100 111 


6 koV 


3.5 (5) GeV 


ELEGANT, [91] 


Nal 


442 m 


10 keV 


5 (8) GeV 


COSME 1, [90] 


Ge 


263 m 


6.5 keV 


5.5 (8) GeV 


CDMS, [112] 


Si 
Ge 


10.6 m 


25 keV 


9.8 (16) GeV 
14 (21) GeV 



simulating the distance a particle travels before scattering, we choose the nucleus 
i it scatters from using the relative probability as in eq. (4.14). 

We take the mass density of the crust to be p — 2.7 g/cm^. To explore the 
sensitivity of the result to the composition of the crust we consider two different 
compositions. First we approximate the crust as being composed of quartz, Si02, 
which is the most common mineral on the Earth and is frequently the primary 
mineral, with > 98% fraction. Then we test the sensitivity of the result by using 
the average composition of the Earth's crust: Oxygen 46.5 %, Silicon 28.9 %, 
Aluminium 8.3 % and Iron 4.8 %, where the percentage is the mass fraction of the 
given element. Our test computer runs showed that both compositions give the 
same result up to the first digit, so we used simpler composition for the computing 
time benefit. Since the DM exclusion window we obtain at the end of this section 



75 



should be very easy to explore in a dedicated experiment, as we show later in the 
text, we do not aim to find precise values of the signal in the underground detector. 

When collisions reduce the DM velocity to less than the Earth's radial escape 
velocity, v^sc — H km/s, DM is captured by the Earth and eventually thermalized. 
CoUisions may also reverse the DM velocity in such a way that the DM particle 
leaves the surface of the Earth with negative velocity: effectively, the DM par- 
ticle gets refiected from the Earth. The majority of light DM particles wind up 
being reflected as is characteristic of diffuse media. The percentage of reflected 
particles proves not to depend on the cross section, as long as the mean free path 
is much smaller than the radius of the earth, but it does depend on DM particle 
mass. Light particles have a higher chance of scattering backward and therefore 
a higher percentage of them are reflected. The initial DM flux on Earth equals 
2.4(1.2) 10^ (1 GeV/mx) cm~^s""-'^, taking standard (conservative) parameter val- 
ues. Table 4.2 shows the fraction of initial flux of DM particles on the Earth which 
are captured and thermalized for various mass values. The fraction is, up to a 
percent difference, independent of whether we make the standard or conservative 
parameter choice. 

For DM particles which are not scattered back to the atmosphere and which 
pass the depth of the detector before faUing below the energy threshold of the 
detector, the scattering in the detector is simulated. For composite targets we 

Table 4.2: The percentage of DM particles incident on Earth which are captured, 
when Xint -C Re- 



mass [GeV] 


2 


4 


6 


10 


100 


thermalized [%] 


21 


30 


36 


46 


94 
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simulate collision with different nuclei with probabilities given as in eq. (4.14). If 
the energy of the recoil is above Ethi we accept the event and record the velocity 
of the particle which deposited the signal. The spectral rate per (kg day keV) is 
then calculated as a sum of rates on the separate elements of the target, as 



where /j is the mass fraction of a given element in the target, px is the local 
DM energy density, /\E is the size of an energy bin of a given experiment and 
< v{a{t)) > is calculated as in (4.16). The signal in the detector falls exponentially 
with (TxN since the energy of DM at a given depth gets degraded as an exponential 
function of the cross section, see [95]. Therefore the hmit on gxn is insensitive to 
small changes in the rate in the detector coming from changes in px', we adopt the 
commonly used value px = 0.3 GeV cm~^ for the local DM energy density. 

We emphasize here that the spectral rate is a function of the relative angle a{t) 
between the direction of the motion of LSR and the position of the detector. This 
angle changes during the day as 



where 6 is the angle between the Earth's North pole and the motion of LSR; a' 
is the angle between the position of the detector and the North pole, and it has 
a value of (90° - geographical latitude). The angle between the LSR motion and 
Earth's North pole is 5 = 42°, so for an experiment around 45° latitude (as for 
Gran Sasso experiment), a' — 45°. Therefore, in the course of a day, the angle 
between the detector and the LSR motion varies in the range approximately 0° to 




(4.18) 



cos a{t) — cos 5 cos a' + sin 5 sin a' sin(a;t + 0o) 



(4.19) 



90°. 



Fig. 4.1 shows the rate R per (kg • day) as a function of time, (4.19), 




i=Al,0 



(4.20) 
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a=400 nb, m^=1.5 GeV. 
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Figure 4.1: The time dependence of the measured rate in underground detectors 
for m—2 GeV and m=1.5 GeV DM candidates. 

calculated for the parameters of the CRESST experiment. We choose 0o so that 
time t=0 corresponds to the moment the detector starts to move away from the 
LSR axis. We see that for these masses the rate is a strong function of the angle of 
the position of the detector with respect to the motion of the LSR, which gives an 
interesting detection signature for detector locations such that this angle changes 
during the day. 

To extract our limits, we average the signal from the simulation dR{t)/dEji 
over one day: 



Since the shape of the spectral rate response is a function of axp in our case 
(because the velocity distribution function at the depth of detector is a function 
of axp due to the interactions in the crust) the extraction of cross section limits is 




(4.21) 
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more complicated than when the rate scales linearly with axp- In the region where 
the window is located, i.e. masses below 2 GeV, we perform the analysis based 
on the fit method used by the CRESST group, [87]. The measured spectrum is 
fit with an empirical function called B. In our case B is the sum of two falling 
exponentials and a constant term, since we expect the signal only in the few lowest 
energy bins. For the fit we use the maximum likelihood method with Poissonian 
statistics in each bin. The maximum likelihood of the best fit, Bq, is Lq. We define 
the background function B' as the difference between the best fit to the measured 
signal, Bq and some hypothesized DM signal S: B' — Bq — S. Following the 
CRESST procedure, we set B' to zero when S exceeds Bq. When (Tq is such that 
the simulated signal S is below the measured signal Bq, B' adds to the signal S, 
completing it to Bq and the likelihood is unchanged. With increasing (Jq, when S 
starts to exceed the function Bq, B' becomes zero, and we calculate the hkelihood 
using S alone in such bins, leading to a new likehhood L. Following the CRESST 
collaboration prescription, ctq excluded at 90% CL is taken to be the value of Uo 
giving InL — InLo = —1.28^/2 [87], since 10% of the area of a normalized Gaussian 
distribution of width a is 1.28cr above than the peak. We show the window obtained 
this way in Fig. 4.4 and for the low mass range, in Figure 4.3. 

For masses higher than 2 GeV we can use a simpler method, since this range of 
masses is already ruled out and our plot is only indicating from which experiments 
the constraints are coming. We calculate the response of the detector for different 
cross section values, and take the limit to be that value for which the simulated 
signal is below the experiment's background. Fig 4.2 shows CRESST background 
together with the simulated response for DM particles with mass nix — 2 GeV 
and 10 GeV and various values of cross section. The limits obtained this way for 
different experiments are given in Figure 4.3. 

The only dark matter detector sensitive to particles with mass ^ 4 GeV is 
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Figure 4.2: The CRESST background and the simulated response of the detector 
for masses rrix — 2 and nix = 10 GeV, and different values of spin independent 

cross sections axp- 
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CRESST. Since it is the only experiment with threshold lower than the threshold 
of the balloon experiment by Rich et al., it extends the existing exclusion window 
for intermediate cross sections. For the CRESST experiment we perform the cal- 
culation using both standard and conservative parameters, because the size of the 
exclusion window is very sensitive to the value of mass threshold, and therefore to 
the parameter choice. For other underground experiments we use only standard 
parameters. In the mass range 5 < m < 14 GeV, the ELEGANT and COSME I 
experiments place the most stringent constraints on a DMIC scenario, since they 
are located in shallow sites a few hundred meters below the ground; see Table 4.1.3. 
Other experiments sensitive in this mass range (e.g. IGEX, COSME II) are lo- 
cated in much deeper laboratories and therefore less suitable for DMIC limits. We 
therefore present limits from ELEGANT and COSME I, for masses 5 to 14 GeV. 
Masses grater than 14 GeV are above the threshold of the CDMS experiment and 
this experiment places the most stringent lower cross section limit due to having 
the smallest amount of shielding, being only 10.6 m under ground. The CDMS 
I had one Si and four Ge detectors operated during a data run. To place their 
limits they used the sophisticated combination of data sets from both types of 
detectors. Due to the large systematic uncertainty on the Si data the Ge data set 
dominates their combined measurements. To be conservative we assume that only 
Ge detectors are present, which reduces the region excluded by CDMS to m > 14 
Gev. Fig 4.3 shows the cross section limits these experiments imply, for masses 
m < 10^ GeV. 

4.1 A Spin-Dependent limits 

In this section we address the case in which DM has a spin dependent interaction 
with ordinary matter. We consider first purely SD interaction and later we consider 
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Figure 4.3: Overview of the exclusion limits for spin independent DM-nucleon 
elastic cross section coming from the direct detection experiments on Earth. The 
limits obtained by using the conservative parameters, as explained in the text, are 
plotted with a solid line; the dotted lines are obtained by using standard parameter 
values and the dashed lines show limits published by corresponding experiments 
or in the case of XQC, by Wandelt et al. [81]. The region labeled with DAMA* 
refers to the hmits pubhshed in [92]. 
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Figure 4.4: The allowed window for ci^^ for a spin independent interaction. The 
region above the upper curve is excluded by the XQC measurements. The region 
below the lower curve is excluded by the underground CRESST experiment. The 
region m > 2.4 GeV is excluded by the experiment of Rich et al. 



83 




Figure 4.5: The allowed Spin Dependent interaction for {Cxp/CxnY — 1- The 
region above the upper curve is excluded by XQC measurements. The region below 
the lower curve is excluded by CRESST. The region m > 2.4 GeV is excluded by 
the balloon experiment of Rich et al. 

the case in which both interaction types are present. We focus on low masses which 
belong to the cross section window allowed by the experiment of Rich et al. 

If the DM has only a spin dependent interaction with ordinary matter, only the 
small fraction of the XQC target with nonzero spin is sensitive to DM detection. 
The nonzero spin nuclei in the target are: Si29 (4.6 % of natural Si), Tei25 (7 %) 
and Hgj^gg, (16.87 %); their spin is due an unpaired neutron. We calculate the spin 
dependent cross section limits from the XQC experiment the same way as for the 
spin independent case, using the new composition of the target. The limiting value 
of the elastic cross section of DM with protons, cr^^ , is shown in Figure 4.5. Since 
the XQC target consists of n-type nuclei, the resulting cross section with protons 
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10 102 

Figure 4.6: usi vs gsd, for CRESST and XQC experiments, for mass mx — 2 
GeV. The region between two curves is the aUowed region. 

is proportional to the {Cxp/CxnY factor as explained in section II. In Figure 4.2 
we use the value {Cxp/CxnY — 1 which is the minimal value this ratio may have. 
We note that the maximal value of the ratio, based on the EMC measurements is 
{Cxp/Cxnf = 500^ and it would shift the XQC limit by a factor 500^ up to higher 
cross sections (substantially extending the allowed window). 

The spin sensitive element in the CRESST target is Al which has an unpaired 
proton in the natural isotope. We assume that the crust consists only of Al, 
since it is the most abundant target nucleus with non-zero spin. In this case the 
model dependence of the C factor ratio drops out in normalizing the underground 
experiment to the proton cross section. 

The window is extended when compared to the purely spin independent DM 
interaction, as shown in Fig. 4.5. This is mostly due to the fact that sensitive part 
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of the XQC target is substantially reduced. 

In Fig. 4.6, for mass rrix — 2 GeV, we plot the asi vs aso limit, assuming 
both types of interactions are present. An interesting feature in the asi vs aso 
dependence is that, when the spin dependent and independent cross sections on 
the target nuclei are of comparable magnitude, screening between two types of 
targets allows cross sections to be higher for the same rate in the detector than in 
the case when only one type of interaction is present. 

4.1.5 Constraint on the fraction of DMIC 

We now turn the argument around and use the XQC data to place a constraint on 
the fraction of allowed DMIC as a function of its elastic cross section. We restrict 
consideration to values of the cross section which are small enough that we do not 
have to treat energy loss in the material surrounding the sensitive components of 
XQC. The maximal fraction DMIC allowed by XQC data p = n%/n^£^M 

can then 

be expressed as a function of cross section, using (4.13) 

P = ^ [A^Si(^'si)(TSi + iVHg((^'Hg)aHg + {vTe)(TTe)] (4.22) 

n-x J -I 

where all quantities are defined as before. 

The mass overburden of XQC can be approximated as [113]: A = 10~^ g/cm^, 
for off- angle from the center of the field of the detector a = (0° to 30°); A = 
10 g/cin, for a = (30° to 100°); and A = 10^ g/cm^ for a > 100°. The center of 
the field of view points toward I — 90°, b — 60° which makes an angle of 32° with 
the detector position direction. Since DM particles are arriving only from above 
the detector, they will traverse either 10 g/cm^ or 10~^ g/cm^ overburden. 

For example, for values of cross section of about 0.7 mb, m = 2 GeV DM 
particles start to interact in the 10 g/cm^ overburden, thus for cross sections above 
this value our simple approach which does not account for the real geometry of the 
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Figure 4.7: The allowed fraction p of DM candidate as a function of DM-nucleon 
cross section. For each mass, p is calculated up to the values of cross sections for 
which the interaction in the mass overburden of the detector becomes important. 



87 



detector, is not applicable anymore. We therefore restrict our analysis to values 
of the cross section for which neglecting the interaction in the overburden is a 
good approximation. In this domain, the allowed fraction of DM falls linearly with 
increasing cross section, as can be seen in equation (4.22) since < vdm > remains 
almost constant and is given by the halo velocity distribution eq. (4.8). The results 
of the simulation are shown in Fig. 4.7, for a spin independent interaction cross 
section. An analysis valid for larger cross sections, which takes into account details 
of the geometry of the XQC detector, is in preparation [114]. 

4.1.6 Future Experiments 

The window for mx ^2.4 GeV in the DMIC cross section range could be explored 
in a dedicated experiment with a detector similar to the one used in the XQC 
experiment and performed on the ground. Here we calculate the spectral rate 
of DM interactions in such detector, in order to illustrate what the shape and 
magnitude of a signal would be. 

In Fig. 4.8 we plot the rate per (kg-day-eV), for a Si detector and DM particle 
masses of mx — 1 and 2 GeV assuming a spin independent interaction. In the 
case of an exclusively spin dependent interaction, the signal would be smaller, 
approximately by a factor f /A^, where / is the fraction of odd-nuclei in the target. 
The calculation is done for a position of a detector for which the signal would be 
the smallest. We assume a short experiment and do not perform averaging over 
time of a day because that would increase the signal. 

The rate scales with cross section; the rate shown in Fig 4.8 is for axp = 2 
/xb, the lower limit on the cross section window from the XQC experiment for 
m — 1 GeV. Since the unshielded muon flux on the ground is of the order of 
2 10^ {w? s)~^ — 2 10^ (cm^ day)~^, an experiment performed on the ground with 
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Figure 4.8: The simulated minimum rate per (kg day eV) calculated with axp — 
2 /ih, for a DM experiment on the ground, versus deposited energy Er in eV, for 
a SI target and for a target with mass number A=100. The sohd hne indicates 
maximal value of the cosmic ray muon background determined based on the total 
muon flux as is used in the text. 
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an array of micro-calorimeter absorbers such as XQC whose target mass is ~ 100 
g, should readily close this window or observe a striking signal. 

4.1.7 Summary 

In § 4.1 we have determined the limits on dark matter in the low mass range 
(m ^ 10 GeV) and with an intermediate cross section on nucleons based on the 
final XQC data and results of underground experiments with low mass threshold. 
We also updated previous limits taking into account newer halo velocity distribu- 
tion. We found that there is an allowed window for DM mass m ^ 2.4 GeV and 
cross section cr ~ //b. Curiously this window overlaps with the mass/cross section 
range expected for the H dibaryon making a possible DM candidate, [40, 66] and 
Chapter 3. We showed that it should be straightforward experimentally to explore 
the window. A signal due to a light DMIC would have strong daily variations 
depending on the detectors position with respect to the LSR motion and therefore 
provide strong signature. 

4.2 The HH Dark Matter- Indirect detection con- 
straints 

B-sequestration scenarios imply the possibility of detectable annihilation of DM 
with anti-baryon number with nucleons in the Earth, Sun or galactic center. The 
rate of H annihilation in an Earth-based detector is the H flux at the detector, 
times cr^^, times (since annihilation is incoherent) the number of target nucleons 
in the detector, 6 x 10^^ per kton. The final products of HN annihilation are mostly 
pions and some kaons, with energies of order 0.1 to 1 GeV. The main background 
in SuperK at these energies comes from atmospheric neutrino interactions whose 
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level is ~ 100 events per kton-yr [115]. Taking <1>^ = -RcapMsii', where Ask i 
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the area of SK experiment and Reap is taken from Table 2.1, the annihilation rate 
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6 10-44cm2_ 

. The total energy release oimH+Bum]^ should give a dramatic signal, so it should 
be possible for SuperK to improve this limit. Note that for the H scenario this 
limit is already uncomfortable, since it is much lower than the effective cross section 
required at freezeout (crg^" = 2.2 10~^^ cm^). However this cannot be regarded 
as fatal, until one can exclude with certainty the possibility that the annihilation 
cross section is very strongly energy dependent. 

Besides direct observation of annihilation with nucleons in a detector, con- 
straints can be placed from indirect effects of H annihilation in concentrations of 
nucleons. We first discuss the photons and neutrinos which are produced by decay 
of annihilation products. The signal is proportional to the number of nucleons di- 
vided by the square of the distance to the source, so Earth is a thousand-fold better 
source for a neutrino signal than is the Sun, all other things being equal. Since 7's 
created by annihilation in the Earth or Sun cannot escape, the galactic center is 
the best source of 7's but do not pursue this here because the constraints above 
imply the signal is several orders of magnitude below present detector capabilities. 

The rate of observable neutrino interactions in SuperK is 



NsK I ^al%^.^dE, (4.24) 



where the sum is over neutrino types, A^sk is the total number of nucleons in 

SuperK, is the spectrum of i-type neutrinos from an H annihilation, a^^j^ is the 
neutrino interaction cross section summed over observable final states (weighted 
by efficiency if computing the rate of observed neutrinos), and is the Ui flux 
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at SK. This last is f^^, the mean effective number of z/j's produced in each H 
annihilation discussed below, times the total rate of H annihilation in the source, 
r^°, divided by 47ri?g, where Rg is the distance from source to SuperK; Rg Re 
for annihilation in Earth. 

In general, computation of the annihilation rate is a complex task because 
it involves solving the transport equation by which DM is captured at the surface, 
migrates toward the core and annihilates, eventually evolving to a steady state 
distribution. However if the characteristic time for a DM particle to annihilate, 
^ann _ ^^aniij^^^^-i^ -g gj^Qp^ Compared to the age of the system, equihbrium be- 
tween annihilation and capture is established (we neglect the evaporation which is 
a good approximation for M^m ~ 0{GeV) and is also more conservative approach) 
so r^'^ equals fcap^g^nR^- Then the neutrino flux, eq. (4.24), is independent of 
"^^jv' because the annihilation rate per H is proportional to it but the equilibrium 
number of ^'s in Earth is inversely proportional to it. For Earth, the equilibrium 
assumption is applicable for a^"" ^ 5 x 10~'*^cm^, while for the Sun it is applicable 
if, roughly, a^^'^^^ 10~^^cm^. For lower annihilation cross sections, transport must 
be treated. 

The final state in HN annihilation is expected to contain A or E and a kaon, 
or S and a pion, and perhaps additional pions. In a dense environment such as the 
core of the Earth, the antihyperon annihilates with a nucleon, producing pions and 
at least one kaon. In a low density environment such as the Sun, the antihyperon 
decay length is typically much shorter than its interaction length. In Earth, pions 
do not contribute significantly to the neutrino flux because 7r°'s decay immediately 
to photons, and the interaction length of 71=*= 's is far smaller than their decay length 
so they eventually convert to vr'^'s through charge exchange reactions; similarly, the 
interaction lengths of K^^s and i^T^'s are much longer than their decay lengths, so 
through charge exchange they essentially all convert to Kg^s before decaying. The 
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branching fraction for production of z/e,^t and Ue,fi from Kg ttI^u is 3.4 x 10""^ 
for each, so /t,. > 2(3.4 x 10~^) for H annihilation in Earth. Since the Sun has a 
paucity of neutrons, any kaons in the annihilation products are typically K'^ and 
furthermore their charge exchange is suppressed by the absence of neutrons. The 
branching fraction for K'^ /^"''^'^t is 63% and the i^^ has 240 MeV if the kaon 
is at rest. If the final states of H annihilation typically contain kaons, then is 
0{1). However if annihilation favors S production, /j, could be as low as ~ 3 - 10~^ 
for production of Pe's and P^^'s above the charged current threshold. Thus the 
predicted neutrino signal from H annihilation in the Sun is far more uncertain 
than in Earth. 

Neutrinos from H annihilation can be detected by SuperK, with a background 
level and sensitivity which depends strongly on neutrino energy and flavor. Taking 
the captured H flux on Earth from Table 2.1, assuming the neutrinos have energy 
in the range 20-80 MeV for which the background rate is at a minimum, and taking 
the effective cross section with which i/'s from the kaon decays make observable 
interactions in SuperK to be 10~^^cm^, eq. (4.24) leads to a predicted rate of excess 
events from annihilations in Earth of r,^sK ~ 2/(kton yr) in the H scenario. This is 
to be compared to the observed event rate in this energy range 3/ (kton yr) [116], 
showing that SuperK is potentially sensitive. If a detailed analysis of SuperK 's 
sensitivity were able to establish that the rate is lower than this prediction, it would 
imply either that the H, H model is wrong or that the annihilation cross section is 
so low that the equilibrium assumption is invalid, i.e., cr^^<2 x 10~^^cm^. The 
analogous calculation for the Sun gives r,^sK ~ 130/,^/ (kton yr) for energies in the 
sub-GeV atmospheric neutrino sample, for which the rate is ~ 35 events/(kton 
yr) [115]^. Thus if f„ were large enough, SuperK could provide evidence for the 
^This estimate disagrees with that of Goldman and Nussinov (GN) [117], independently of the 
question of the value of . GN use an H flux in the solar system which is eight times larger than 
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H, H scenario via energetic solar neutrinos, but the absence of a solar neutrino 
signal cannot be taken as excluding the H scenario, given the possibility that 

/. < 10-=^. 

Fortunately, there is a clean way to see that the DM cannot contain a sufficient 
density of ^'s to account for the BAU. When an H annihilates, an energy rriH + 
Bhttln is released, almost all of which is converted to heat. Uranus provides a 
remarkable system for constraining such a possibility, because of its large size 
and extremely low level of heat production, 42 ± 47 erg cm~^s~^, (Uranus internal 
heat production is atypically small, only about a tenth of the similar sized planet 
Neptune), [118]. When annihilation is in equilibrium with capture as discussed 
above, the power supplied by annihilation is P^""^ = f^^^(^xij^x + BxTTIn). For 
the H, fti 0.2 as for Earth, so the heat flux generated in Uranus should be 
470 erg cm~^s~^, which definitively excludes the H, H scenario. 

4.3 Conclusion 

In this section we have shown that the H di-baryon could evade experimental 

searches if it is compact and tightly bound. It would not bind to nuclei and 

therefore the anomalous mass isotope experiments would not be sensitive to its 

existence. For masses m^f <2.05 GeV it would also be cosmologically stable. As 

such it could potentially offer an explanation of DM problem within the Standard 

Model. We showed that the H alone could not be produced in the DM abundance 
our value in Table 2.1 from integrating the normal component of the halo velocity distribution, 
due to poor approximations and taking a factor-of-two larger value for the local DM density. 
We include a factor 0.35 for the loss of I'^'s due to oscillation, we account for the fact that only 
neutrons in SuperK are targets for CC events, and we avoid order-of-magnitude roundup. Note 
that the discussion of the particle physics of the H in [117] applies to the case of an absolutely 
stable H, which we discussed but discarded in [66]. 
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through thermal decouphng from SM particles. In the B-sequestration scenarios, 
HH could be produced with proper abundance in the early universe at a decou- 
phng temperatures of around 85 MeV. We find that mass and cross section ranges 
expected for the H is not ruled out by current DM direct detection experiments, 
but that the H, H scenario could be ruled out through the heat production due to 
HN annihilation in Uranus. 
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Chapter 5 

New Particle X with Baryon 
Number 

5.1 Particle properties of X 

We now turn to the possibility of a new, light fundamental particle with Bx — 1 
and mx ^4.5 GeV. Such a low mass suggests it is a fermion whose mass is pro- 
tected by a chiral symmetry. Various dimension-6 interactions with quarks could 
appear in the low energy effective theory after high scale interactions, e.g., those 
responsible for the family structure of the Standard Model, have been integrated 
out. These include 

K{Xh d'^c - Xc M) + h.c, (5.1) 

where the b and c fields are left-handed SU(2) doublets, combined to form an SU(2) 
singlet, d'^ is the charge conjugate of the SU(2) singlet field du, and k — fl'^/A^, 
where A is an energy scale of new physics. The suppressed color and spin indices 
are those of the antisymmetric operator O" given in equation (10) of ref. [119]. 
The hypercharge of the left-handed quarks is +1/3 and that of dR is -2/3, so the 
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X is a singlet under all Standard Model interactions and its only interaction with 
fields of the Standard Model are through operators such as eq. (5.1). Dimension-6 
operators involving only third generation quarks can be constructed; supplemented 
by W exchange or penguins, they could also be relevant. Note that k, is in general 
temperature dependent and we denote its value today (at freezeout) by Kq and Kfo 
respectively. 

Prior to freezeout, X's stay in equilibrium through scattering reactions like 

d + X^bc. (5.2) 

The coupling k, in eq. (5.1) is in general complex and a variety of diagrams involving 
all three generations and including both W exchange and penguins contribute to 
generating the effective interaction in eq. (5.2), so the conditions necessary for a 
sizable CP-violating asymmetry between and (7^° are in place. 
An interaction such as eq. (5.1) gives rise to 

_ 1 nixmhrricTfo 

For the freezeout of X to occur at the correct temperature (see Table 2.1), Kfo 
10~^GeV~^ is needed. This suggests an energy scale for new physics of A<10 
TeV, taking dimensionless couplings to be < 1. 

X particles can decay via production of bed quarks, but this state is too massive 
to be kinematically allowed. For an X particle mass of 4.5 GeV the decay is off- 
shell, with a W exchange between b and c quarks, giving: 

X csd (5.3) 

The matrix element for this transition is proportional to 

M ^ «:^^|HcKs| / d'k^l^^ (5-4) 
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The integral over the loop momentum gives \n{h./Myj) and the diagram is loga- 
rithmically divergent The decay rate of X today can be estimated as: 



where gw is the electroweak SU (2) gauge coupling. The condition tx ^ Tuni^erse 
places a constraint on the value of the X coupling today, kq < 10"^'^ GeV~^. Thus 
for X to be a valid dark matter candidate, its coupling to ordinary matter needs to 
have a strong temperature dependence changing from 10~^ GeV~^ at a temperature 
of ~ 200 MeV, to 10~^°GeV~^ or effectively zero at zero temperature. If the 
interaction in eq. (5.1) is mediated by a scalar field r] with couplings of order 
1, its effective mass m,, should vary from 10 TeV at 200 MeV to 10^° TeV at 
zero temperature. The most attractive way to do this would be if it were related 
somehow to a sphaleron-type phenomenon which was allowed above the QCD or 
chiral phase transition, but strongly suppressed at low temperature. We do not 
attempt that here and instead display two "toy" examples of models where the 
desired dramatic change in k occurs. Let the dominant contribution to the t] mass 
be due to the VEV of another scalar field a which has baryon and other quantum 
numbers zero. The VEV of a can be rapidly driven from zero to some fixed value 
resulting in the desired mass change in 77 by several possible mechanisms. The 
simplest requires only one additional field with the zero temperature interaction 



r ~ m\i<lg^\VbcVcs 



(5.5) 



The global minimum of this potential at zero temperature is at 




(5.6) 



The mass of the field 77 in this scenario equals 



^2^3^ = ^/{a^/ai)m^ ~ 10^° TeV. 



(5.7) 
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At higher temperature, one loop corrections contribute to the potential and intro- 
duce a temperature dependence [120, 121]: 

Vioop = g T a H T 7] . (5.8) 

The new condition for the minimum of the potential in the a direction becomes 

^ m;-r'(2a, + a3)/3_ 
4q;i 

and for temperatures higher than the critical value 

Zai + as 

the potential has a unique minimum at (77) = 0, (a) — 0. Condition (5.7) together 
with TcR ~ 200 MeV imphes the relation 

as ~ 10^ (5.11) 

This large difference in coupling looks unnatural. Fine tunning can be avoided at 
the price of introducing a second helping scalar field as in Linde's hybrid inflation 
models [122]: 

V{a,<t>) = i-(M^-AaY + ^ + ^0V (5.12) 

M^_fM^_g^\ , \ 

4A V2 27^ 4^ 2 

The potential is such that, for values of > M/g, its minimum in the a direction 

is at (a) = 0. The evolution of fields (j) and a as the universe expands would be as 
follows. At high temperatures we assume that the field is at a high value, 

, M 

<t>>—- (5-13) 
The equation of motion of the field in an expanding Universe is 

+ 3/^0 + \/'(0) = 0. (5.14) 
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The solution for radiation dominated universe, where Hubble constant scales as 
H= l/{2t) is 

^ Ci + C'^{mt)-^/^, mt^O (5.15) 

where J and Y are spherical Bessel functions, and (p and oscillates for sufficiently 
large mt. As rolls down the potential it becomes equal to (f) — ^ the 
symmetry breaking occurs. The potential develops a minimum in the a direction 
and the value of a field tracking the minimum becomes 

M = - — ^^-^^^ 

As the temperature drops further goes to zero on a time scale 1/m and the VEV 
of (7 goes to its asymptotic value 

(a) = ^. (5.17) 

From the condition on the value of coupling of X today, with k,q ~ 1/"^,,^ ~ l/{a)^ 
and assuming A ~ 1 we get the value of parameter M of M > 10^° GeV. We can 
place a constraint on ^ (0) from a condition that the energy density in the field 0, 
V{(f)) — m^0^/2 should be less than the energy density in radiation at temperatures 
above ~ 200 McV. Since the field (p rolls down slowly as t^^^^ and ~ while 
radiation scales down more rapidly, as prad ~ ~ it is enough to place the 
condition on the value of the energy density in the field cf) at 200 MeV, Prad^ Pcf,- 
This leads to the condition that the mass parameter for cf) must satisfy m < 10"^ 
eV. Requiring that mi < 1 at T = 200 MeV, to avoid the oscillation phase, sets the 
condition m < 10^^^ eV. Achieving such low masses naturally is a model-building 
challenge. Work on the cosmological implications of these scenarios is in progress 
and will be presented elsewhere. 
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5.2 Constraints on XX DM 



The approach presented here to solve the DM and BAU puzzles at the same time, 
with baryonic and anti-baryonic dark matter, can run afoul of observations in 
several ways which we now check: 

5.2.1 Direct detection constraints 

The scattering cross sections (Txn and (Txn must be small enough to be compat- 
ible with DM searches - specifically, for a 4 GeV WIMP. If Xs interact weakly 
with nucleons, standard WIMP searches constrain the low energy scattering cross 
section aoM = (c"^jv ^^xn)/^^ + Table 2.1 gives the capturing rate of X by 
the Earth, Reap- The capturing rate is obtained using code by Edsjo et al. [29] 
which calculates the velocity distribution of weakly interacting dark matter at the 
Earth taking into account gravitational diffusion by Sun, Jupiter and Venus. It is 
not possible to use the upper limit on kq from requiring the X lifetime to be long 
compared to the age of the Universe, to obtain o"^jv{xjv} 'without understanding 
how the interaction of eq. (5.1) is generated, since it is not renormalizable. A naive 



guess 



m'irm'i 



'^xNixm - ^ A ^^^^^^y (5.18) 

is well below the present limit of ~ 10~^*cm^ for a 4 GeV particle, even using 
the maximum allowed value of kq, but the actual value depends on the high-scale 
physics and could be significantly larger or smaller. 

5.2.2 Indirect constraints 

A second requirement is that the annihilation of X with matter does not produce 
observable effects. If eq. (5.1) is the only coupling of X to quarks and ~ 10~^° 
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GeV~^, the effects of anniliilation in Eartli, Sun, Uranus and tlie galactic center 
are unobservably small. In this case the very stability of the X implies that its 
interaction with light quarks is so weak that its annihilation rate in the T — Q 
Universe is essentially infinitesimally small. The cross section for the dominant 
annihilation processes is governed by the same Feynman diagrams that govern X 
decay so that dimensional arguments lead to the order of magnitude relation 

^Fjv ~ ^x^^x^ - 10-^^(30Gyr/Tx) cm^ (5.19) 

For completeness, in the rest of this section we will discuss the indirect limits 
which could be placed from annihilation experiments on the X4 DM, although, as 
we have seen, the expected X^ cross sections are much smaller than what current 
limits could demand. 



Direct detection of annihilation. The discussion in this subsection is similar 
to the analysis of H as DM in B-sequestration models in §4.2. The rate of X 
annihilation in SupcrK detector is, analogously to the case of H (see §4.2), the X 
flux at the detector times c^^, times the number of target nucleons in the detector 
(since annihilation is incoherent) and has the value 



a 



X 



10-45 cni2 



(kton yr) ^ . 



(5.20) 



The X signal is lower than the background if (J'xno ^ 2 x 10 "^^cm^, which is 
readily satisfied as we have seen above. 



Indirect detection of annihilation: Besides direct observation of annihilation 
with nucleons in a detector, constraints can be placed from indirect effects of X 
annihilation in concentrations of nucleons. 
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Neutrinos from X annihilation can be detected by SuperK, with a background 
level and sensitivity which depends strongly on neutrino energy and flavor. The 
rate of observable neutrino interactions in SuperK is given by eq. (4.24). We 
will distinguish two contributions to the DM annihilation rate resulting in 
a neutrino signal, 1) annihilation of a total DM flux occurring while DM 

passes through the Earth and 2) annihilation rate of the small percentage of DM 
particles that are gravitationally captured due to scatter from nuclei in the Earth 
and therefore eventually settle in the Earth's core. In the first case the annihilation 
rate can be estimated as rj^^ ~ Ne $^■^(7^'^ where Ne is the number of nucleons 
in the Earth. We assume that the annihilations are spread uniformly in the Earth 
and that /j^ neutrinos is produced in each annihilation. We calculate the signal in 
SuperK due to neutrino interaction in the detector, taking the effective cross section 
with which i/'s from the kaon decays make observable interactions in SuperK to 
be 10~^^cm^, in eq. (4.24) and we get the signal in SK of 



^ann 
" X 



10-45 cm- 



(kton yr)-^ (5.21) 



The annihilation of DM in the Earth and subsequent neutrino production therefore 
do not produce detectable signal for fi, < 10^. 

In general, computation of the annihilation rate in the case of captured DM is 
a complex task because it involves solving the transport equation by which DM is 
captured at the surface, migrates toward the core and annihilates, eventually evolv- 
ing to a steady state distribution. However in equilibrium, and when neglecting the 
evaporation, capturing rate equals to the annihilation rate, F^"^ = /cap$x47ri?|;, 
see § 4.2. Then the neutrino flux in eq. (4.24) is independent of cr^jlj^, as we have 
seen in § 4.2. 

Taking the captured X flux on Earth from Table 2.1, eq. (4.24) leads to a 
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predicted rate of the neutrino interaction in SK of 



X 



10-45 cm2 



(yrkton) ^. 



(5.22) 



The analogous calculation for the Sun gives even smaller rates for energies in the 
sub-GeV atmospheric neutrino sample. The most stringent constraint in this B- 
sequestration model therefore comes from DM annihilation in SuperK, eq. (5.20), 
and it is safe for cross sections of interest. 
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Chapter 6 

Summary and Outlook 



In this thesis we addressed several problems related to the nature of dark matter: 
we summarize here briefly our results. 

The existence and properties of the H dibaryon. Since it was predicted 
to be a strong interaction stable 6 quark state, this particle raised a huge interest 
because it pointed to possibility of new type of hadrons. As we have summarized 
in §3, extensive experimental effort has been made in trying to produce and detect 
it. Recently, experiments on double A hypernuclei claimed to rule it out. In the 
work presented in the thesis, we show that, if suffitiently compact {th^^/'^tn) 
the H formation time from a double A system could be longer than the single A 
decay lifetime and therefore the double A experiments would be insensitive to its 
existence. Furthermore, we discover that, for an even smaller H (rj^ < l/3rAr), 
with a mass smaller than rnH'^TriN + m\, the H could be cosmologically stable. 
We also find that, for the reasonable range of values of coupling to a meson and 
glueball the H would not bind to nuclei, and therefore the anomalous mass isotope 
experiments cannot rule out its existence. 
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Can the if be a DM candidate? Is a Standard Model Dark Matter 
ruled out? Given that the if is a neutral particle which could be sufficiently 
long lived, and having the virtue that it was predicted within the SM, we posed the 
question weather it could be DM. We analyzed the results from DM experiments 
sensitive in the H expected mass and cross section ranges in §4. This region 
of the parameter space was not analyzed before since only the new generation of 
underground experiments reached such low mass region. We also analyzed the final 
data from the X-ray experiment XQC, which proved to exclude a large part of the 
parameter space for the low masses and high cross sections. Surprisingly, there is an 
allowed window in DM exclusion region. We show that the window should be easy 
to close with a dedicated experiment. Given that the H would not bind to nuclei 
and that it would be nonrelativistic after the QGD phase transition, the H would 
be a cold DM candidate allowed by experiments. The production mechanism in the 
early Universe turns out to be problematic, since the H could not be produced in 
sufficient amount by the standard mechanism of thermal production. The reason 
for this is that in order to be abundant enough it would need to stay in the 
equilibrium until low temperatures (T ~ 15 MeV), when all the strange particles 
needed for its production have already decayed. 

Can Dcirk Matter Ccirry baryon number and be the answer to the baryon 
asymmetry of the Universe? We worked on a scenario in which dark matter 

carries (anti)baryon number and offers a solution to the baryon asymmetry prob- 
lem. We analyzed two concrete models, the if, H DM and the new Beyond the 
Standard Model particle X as candidates for the model. For the if, H scenario we 
already checked that DM detection experiments allow for its existence and that it 
has the correct particle properties to be undetected and long lived. In this scenario 
the new set of constraints with respect to the H DM comes from the annihilation 
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of in regions with high concentration of nucleons. The H can successfully evade 
constraints from the direct detection of annihilation in SuperK, and the detection 
of neutrinos produced by its annihilation in the Sun and the Earth. However, 
the heat production in Uranus, which is a planet with an anomalously low internal 
heat production, is lower than the heat that would be produced by the annihilation 
of captured Hs. This excludes a HH dark matter. The other scenario, involving 
the new particle X turns out to be safe from the above constraints. Its stability 
requires that its coupling to quarks should have a temperature dependence, and we 
analyze two models which could provide the change in the coupling. It also follows 
that the value of the coupling today is such that the X is virtually undetectable 
by current experiment. 
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Appendix A 

MC simulation for DM heavier 
than 10 GeV 



We assume the following function for the form factor, as explained in Section 4.1.1, 

F\q^) ^ exp-ToiiR)\ (A.l) 

where q is momentum transfer and R is the nuclear radius. For a particle mov- 
ing with a given velocity v, the mean free path to the next collision is obtained 
using the cross section atot which corresponds to a{q) integrated over the avail- 
able momentum transfer range, from zero to qmax, where q^ax — '^''^NER^rnax ^-nd 
ER,max = 'lli^lmAivjcf: 

_ \l-F\q-)dq- 

After a particle travels the distance calculated from the mean free path described 
above, the collision is simulated. The momentum transfer of a collision is deter- 
mined based on the distribution given by the form factor function, as in the usual 
Monte Carlo method procedure 

^~ j^-^F^{q^)dq^' ^ "^^ 
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where p is a uniformly distributed random number from to 1. Once the momen- 
tum transfer of the coUision is determined, the recoil energy of the nucleus, Er, 
and the scattering angle of the collision, 6cm i are uniquely determined. 

We repeat this procedure while following the propagation of a particle to the 
detector. If the particle reaches the detector we simulate the coUision with target 
nuclei. For each collision in the target, the energy deposited in the detector Er is 
determined as above. For each particle i the energy transfer determines the cross 
section with target nuclei as crxAi{ER) — crxA,oF"^{ER). The rate in the detector 
is found as in equation (4.18) with the only difference that in this case the sum 
runs over < v{a{t))axAi >i instead of depending only on v{a{t)). 
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Appendix B 

Relative probability for scattering 
from different types of nuclei 



The probability P{x+dx) that a particle will not scatter when propagating through 
a distance x+dx, equals the probability P{x) that it does not scatter in the distance 
X, times the probability that it does not scatter from any type i of target nuclei in 
the layer dx: 




By solving this differential equation one gets the probability that a particle will 
travel a distance x in a given medium, without scattering, 



The probability for scattering once and from a given nuclear species i in the 
layer {x, x + dx), is proportional to the product of probabilities that a particle will 
not scatter in distance x and that it will scatter from species of type i in dx: 



P{x) = e 



(B.2) 



fi{x)dx = e 




(B.3) 
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The probability that a particle scatters once from any species in a dx layer is 
the sum of the single particle probabilities ^ fi{x)dx, where 

/"OC 

/ y2ux)dx = i. (B.4) 

In the simulation wc want to generate the spectrum of distances a particle trav- 
els before scattering once from any of elements, using a set of uniformly distributed 
random numbers. We can achieve this by equating the differential probability for 
scattering to that of a uniformly distributed random number, 

J2Mx)dx^dR (B.5) 

After integrating 

/ Y] M^) dx= dR (B.6) 

Jo Jo 

we get for the distribution of scattering distances x 

X = -Xeff InR (B.7) 
The relative frequency of scattering from a nucleus of type i, is then given by 



/ fi{x)dx = = — B.8 
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